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INTRODUCTION 

The  purpose  of  the  work  funded  by  this  grant  is  to  develop  a  new  type  of  dedicated  system  for 
diagnostic  breast  imaging  that  can  simultaneously  obtain  structural  information  from  a  digital 
mammography  detector  and  functional  information  from  a  high  resolution  gamma  imaging  detector. 

The  overall  goals  of  the  project,  as  stated  in  the  original  application,  are: 

a)  To  upgrade  and  optimize  our  test  dual  system, 

b)  To  compare  its  performance  over  a  broad  range  of  phantom  imaging  tasks  to  that  of  the 
combination  of  film-screen  mammography  and  scintimammography,  and 

c)  To  perform  a  limited  clinical  study  in  which  the  results  of  the  dual  modality  study  are 
compared  to  the  results  of  biopsy. 

BODY 

As  taken  from  the  original  Statement  of  Work,  the  tasks  scheduled  to  begin  and/or  be  completed 
during  the  project  period  are  as  follows: 

Task  1.  Upgrade  digital  mammography  and  gamma  detector  components  of  test  system.  Months  1-6 

a)  Replace  current  6  cm  x  6  cm  spot  detector  with  20  cm  x  30  cm  full  field  digital 
mammography  (FFDM)  detector 

b)  Characterize  imaging  properties  of  FFDM  detector,  including  linearity,  dynamic  range, 
modulation  transfer  function  (MTF),  noise  power  spectrum  (NFS),  noise  equivalent  quanta 
(NEQ),  detective  quantum  efficiency  (DQEX  geometric  linearity,  and  artifacts. 

c)  Replace  6  cm  x  6  cm,  1.25  mm  pixel  detector  with  8  cm  x  8  cm,  2.25  pixel  detector 

d)  Characterize  detector  imaging  properties  including  sensitivity,  line  spread  function  (LSF), 
energy  resolution,  geometric  linearity,  and  uniformity  of  response. 

e)  Modify  co-registration  software  to  accommodate  new  detector  parameters 
Task  2.  Collimator  optimization.  Months  3-6 

a)  Calculate  collimator  efficiency,  resolution  for  the  range  of  object  distances  relevant  in  the 
dual  system  for  several  potential  materials  (e.g.  tantalum,  lead  or  tungsten). 

b)  Fabricate  and  evaluate  small  test  collimators  covering  a  limited  fraction  of  the  FOV, 
including  slant  hole  collimator. 

c)  Build  optimized  high  resolution  and  optimized  high  efficiency  collimators  for  dual 
system. 

Task  3.  Phantom  development.  Months  7-12 

a)  Calculate  useful  range  of  phantom  parameters  (thickness,  lesion  sizes,  etc)  based  on  x-ray 
and  gamma  ray  attenuation,  detector  resolution,  and  scatter  measurements. 

b)  Fabricate  fillable  contrast-detail  phantoms,  chamber  phantoms  with  removable  simulated 
lesions,  fillable  vessels  of  various  thicknesses  to  provide  background  radioactivity,  etc. 

c)  Purchase  breast  equivalent  mammography  block  phantoms  with  imbedded  simulated 
lesions  (CIRS,  Newport  News,  VA) 

d)  Evaluate  dual  system  performance  using  permutations  of  stacked  phantoms,  scatter  free 
conditions 
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e)  Repeat  the  measurements  of  d),  but  incorporating  the  anthropomorphic  torso  phantom  to 
create  realistic  scatter  conditions. 

Taj'fc  Quantification  of  the  effects  of  compression,  Months  7-18 

a)  Use  the  phantoms  developed  in  Task  3  to  measure  SNR,  target-to-background  ratio  while 
varying  the  thickness  of  interposed  phantom  material  and  total  phantom  thickness 

Task  5.  Comparison  of  dual  system  performance  with  diagnostic  mammography  and 
scintimammography.  Months  19-30 

a)  Repeat  the  measurements  of  Tasks  3d-e  using  a  high  resolution  clinical  Anger  camera 
with  10”  field  of  view  (Technicare),  and  clinical  screen-film  mammography  systems  (Lor ad 
M-IV,  Siemens  Mammomat). 

b)  Quantify  lesion  contrast  in  gamma  and  x-ray  images  (using  a  film  densitometer),  for  both 
dual  and  clinical  systems. 

c)  Perform  reader  study  of  contrast-detail  images  from  dual  system  and  from  clinical 
systems.  Readers  will  be  physicists,  technologists,  and  radiologists  at  UVa.  Digital  x-ray 
images  will  be  printed  to  laser  film  for  the  study. 

Task  6.  Evaluation  of  multiple  gamma  and  x-ray  views.  Months  24-30 

a)  Modify  phantoms  by  including  blocks  (modules)  that  are  composed  of  heterogeneous 
materials  at  several  layers  to  simulate  superimposed  beast  structure 

b)  Modify  breast  support  and  compression  paddle  to  be  fixed  and  independent  of  the 
detector  gantry 

c)  Use  multiple  views,  single  phantom  orientation.  For  each  view,  obtain  corresponding  x- 
ray  and  gamma  ray  images. 

d)  Evaluate  the  efficacy  of  3-d  localization  and  removal  of  superimposed  structures  via 
tomosynthesis,  3-d  reconstruction. 

Task  7.  Patient  study  on  dual  modality  system.  Months  24-36 

a)  Update  IRB  to  reflect  protocol  found  most  advantageous  (i.e.  number  of  views  per 
compression)  for  lesion  conspicuity 

b)  Measure  uncompressed  (pendulant)  breast  thickness,  compressed  thickness,  and 
compressive  force  for  each  study 

c)  Conduct  reader  study  of  dual  images  with  two  board-certified  radiologists  (UVa). 
Correlate  reading  results  with  pathology  report. 

Task  8.  Publication  of  results.  Months  34-36 
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Accomplishments 

As  detailed  in  our  first  year  annual  report,  Task  1  was  completed  and  resulted  in  a  publication  in 
Medical  Physics  (Williams  et  al.  1999).  The  results  of  Task  2  have  also  been  published,  as  part  of  a 
publication  in  TF.F.F.  Transactions  on  Nuclear  Science  (Majewski  et  al,  2001).  Task  3,  the 
development  of  phantoms,  was  also  essential  for  the  results  published  in  the  latter  manuscript,  as 
were  studies  done  under  Task  4  investigating  the  effect  of  breast  compression  on  lesion  detectability 
using  phantoms.  As  described  in  our  report  of  July  31, 2002,  we  have  completed  analysis  of  the 
dependence  on  compressed  breast  thickness  on  the  scatter  to  primary  ratio,  using  data  obtained  in 
our  clinical  study  (Task  7).  These  results  were  presented  at  the  2002  IEEE  Nuclear  Science 
Symposium  and  Medical  Imaging  Conference  in  Newport  News,  VA,  Nov.  10-16, 2002,  and  have 
been  published  in  the  IEEE  Transactions  on  Nuclear  Science  (Williams  et  al.,  2003).  The  results 
show  the  measured  scatter  fraction  (fraction  of  total  counts  in  a  given  region  of  interest  that  are  due 
to  scattered  gamma  rays)  as  a  function  of  compressed  thickness.  The  data  show  that  compression  is 
likely  to  reduce  the  amount  of  detected  scatter. 

As  we  described  in  our  July  31, 2002  report,  we  have  decided  that  Task  5,  as  originally  proposed, 
will  not  be  undertaken  because  in  the  years  since  that  task  was  proposed,  it  has  become  more  clear 
that  it  was  inappropriately  suggested  since  the  scanner  is  being  used  in  lesion  characterization  as 
opposed  to  lesion  detection.  Lesion  characterization  requires  accurate  assessment  of  the  relative 
radioactivity  concentration  in  the  lesion  and  surrounding  healthy  tissue.  Quantitation  of  lesion 
radiotracer  uptake  requires  corrections  for  the  effects  of  camera  spatial  resolution,  gamma  ray 
attenuation,  and  scatter.  The  first  two  of  these  phenomena  are  functions  of  the  lesion-to-collimator 
separation,  and  thus  correction  requires  knowledge  of  the  lesion  depth.  Therefore,  we  have  invested 
effort  in  a)  multiple  view  imaging  (originally  seheduled  as  Task  6),  and  b)  development  of 
quantitation  tools. 

In  regard  to  task  a),  we  described  in  our  July  2002  report  the  hardware  modifications  that  have  been 
made  to  a  second  upright  mammography  unit  to  permit  multiple  images  of  the  breast  to  be  obtained 
with  fixed  breast  compression.  This  setup  permits  highly  reproducible,  computer  controlled 
positioning  of  the  x-ray  tube  and  a  small  field  of  view  (5  cm  x  5  cm)  digital  x-ray  detector.  We 
have  also  mounted  a  larger  x-ray  detector  (20  cm  x  30  cm  field  of  view),  on  the  motorized  gantry 
arm.  An  abstract  describing  early  results  of  multiple-view,  dual  modality  imaging  using  phantoms 
has  been  accepted  for  presentation  at  the  IEEE  Nuclear  Science  Symposium  and  Medical  Imaging 
Conference,  October  16-22, 2004,  Rome,  Italy. 

Regarding  task  b)  we  have  developed  a  mathematical  model  whose  purpose  is  to  calculate  changes 
in  lesion  image  contrast  as  a  function  of  the  location  within  the  breast  of  the  lesion,  breast  thickness, 
and  gamma  camera  position.  The  results  of  this  work  have  been  presented  at  several  international 
meetings  (the  First  International  Symposium  for  Functional  Breast  Imaging  in  Rome,  Italy  R'dine, 
Italy,  April  18-21, 2001.,  the  IEEE  Nuclear  Science  Symposium  and  Medical  Imaging  Conference, 
Norfolk,  VA,  November  12-16,  2002,  and  the  7“*  International  Workshop  on  Digital 
Mammography,  June  18-21, 2004,  Durham,  NC.),  and  has  led  to  three  publications  (see  below). 

We  are  completing  the  reader  study  associated  with  our  clinical  trial.  Altogether  we  have  imaged  26 
volunteers. 

Summary 

To  summarize,  we  have  now  completed  all  of  the  Tasks  originally  proposed  except  Task  5,  which 
we  have  learned  is  inappropriate,  and  part  c)  of  Task  7,  the  reader  study.  We  have  nearly  finished 
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♦ 

the  reader  study,  with  one  reader  remaining.  This  project  has  resulted  thus  far  in  nine  publications 
and  eleven  scientific  presentations.  The  completion  of  the  reader  study  will  result  in  at  least  one 
more  of  each.  Some  of  the  recent  articles  are  included  in  the  Appendix. 


KEY  RESEARCH  ACCOMPLISHMENTS 

•  Upgraded  dual  modality  x-ray  and  gamma  ray  detectors,  with  measured  performance  meeting 
design  specifications,  and  upgraded  the  dual  modality  software 

•  Optimized  gamma  camera  collimator 

•  Developed  new  types  of  phantoms  for  evaluation  of  factors  determining  lesion  contrast  in  breast 
scintigraphy 

•  Quantified  the  effects  of  breast  compression  on  lesion  contrast  in  breast  scintigraphy 

•  We  have  conducted  studies  leading  to  the  use  of  an  improved  gamma  camera  scintillator 
material  permitting  better  scatter  rejection 

•  We  have  evaluated  the  effects  of  breast  compression  on  lesion  contrast  in  breast  scintigraphy 

•  We  have  developed  a  mathematical  model  that  relates  gamma  camera  resolution  and  lesion  size 
and  depth  to  image  contrast 

•  We  have  begun  modifications  of  a  second  mammography  gantry  to  permit  multiple  views  of  the 
breast  with  both  imaging  modalities 

•  We  have  investigated  methods  for  scatter  correction  in  the  scintimammography  images,  and 
have  characterized  the  relationships  between  the  amount  of  scatter  and  the  location  within  the 
breast  and  overall  breast  thickness. 

•  We  have  obtained  very  promising  preliminary  results  in  our  clinical  evaluation  of  the  dual 
modality  scanner,  and  are  now  completing  a  reader  study  to  assess  the  added  value  of  using  dual 
modality  scanning  to  aid  in  diagnostic  workup. 

•iv/n'- 
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CONCLUSIONS 

We  have  developed  a  new  type  of  dual  modality  system  for  breast  cancer  diagnosis.  The  device 
constitutes  an  integrated  system  for  combining  functional  and  structural  information  regarding 
normal  and  pathological  breast  structures.  The  system  has  demonstrated  the  ability  to  provide 
unambiguous  correlation  between  structures  in  the  stmctural  (x-ray)  and  functional  (gamma  ray) 
images.  Upgrading  of  scanner  components  and  associated  software  has  resulted  in  a  more  accurate 
and  versatile  scanner.  Basic  studies  using  phantoms  have  pointed  to  the  value  of  mild  breast 
compression  for  increasing  lesion  contrast.  In  addition,  our  studies  have  pointed  to  the  importance 
of  optimization  of  image  acquisition  strategy  based  on  the  location  of  the  lesion  within  the  breast. 
Characterization  of  the  scatter  radiation  contributing  to  the  gamma  image  has  shown  a  distinct 
correlation  between  breast  thickness  and  the  amount  of  scatter.  Early  results  of  a  clinical  evaluation 
of  the  scanner,  using  histologic  findings  as  ground  truth,  are  very  promising.  Ongoing  work, 
including  multiple  view  acquisition,  is  likely  to  further  enhance  the  diagnostic  capabilities  of  the 
system.  It  is  anticipated  that,  compared  to  current  non-integrated  imaging  procedures,  breast  lesions 
that  are  smaller  and  of  a  wider  variety  can  ultimately  be  reliably  characterized. 
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Phantom  Study  of  Radiotracer  Concentration 
Quantification  in  Breast  Scintigraphy 

Mark  B.  Williams,  Member,  IEEE,  Mitali  J.  More,  Deepa  Narayanan,  Stan  Majewski,  Ben  Welch,  Randolph  Wojcik, 

and  Douglas  A.  Kieper. 


Ahsii^act — The  autliors  are  developing  a  breast  imaging  system 
that  combines  digitai  X-ray  mammography  and  gamma  emission 
scintigraphy  in  an  integrated  unit.  One  potential  value  of  such 
a  combinatioh  is  the  ability  to  quickly  obtain  information  about 
breast  shape  and  lesion  position  within  the  breast  via  X-ray  and 
then  to  use  this  information  to  perfohn  corrections  to  the  gamma 
image(s).  In  particular,  the  effects  of  gamma  ray  attenuation, 
gamma  camera  resolution,  and  scatter  can  be  taken  into  account 
so  as  to  better  estimate  the  difference  between  the  image  contrast 
and  true  radioactivity  contrast.  The  authors  present  results  of  a 
phantom  study  designed  to  test  how  well  the  lesion-to-background 
contrast  can  be  predicted  by  this  method.  The  preliminary  results 
indicate  that  the  effects  of  attenuation,  partial  volume  averaging, 
and  scatter  on  lesion  image  contrast  can  be  modeled  and  suggest 
that  the  model  might  be  usefully  employed  in  the  dual  modality 
breast  scanner  to  correct  the  breast  scintigram  for  these  effects. 

Index  Terms — Digital  mammography,  multimodal  imaging, 
radioactivity  quuutiJicullun,  selntimammography,  stereotactic 
imaging. 


L  Introduction 

Breast  scintigraphy  (scintimammography)  is  being  eval¬ 
uated  as  an  adjunct  modality  for  X-ray  mammography, 
because  it  provides  functional  information  complimentary  to 
the  structural  information  found  in  the  mammograpbic  study. 
We  are  developing,  a  system  for  dual  modality  breast  imaging, 
.in  which  digital  mammography  images  and  breast  scintigraphy 
images  are  obtained  in  quick  succession,,  with  the  breast  held 
under  mild  compression.  The  X-ray  and  gamma  ray  detectors 
are  mounted  on  a  common  gantry,  and  the  fact  that  the  breast 
remains  stationary  during  X-ray  and  gamma  ray  image  acquisi¬ 
tion  permits  unambiguous  correlation  between  the  lesion  loca¬ 
tion  in  the  two  images. 

For  a  given  true  differential  radiotracer  concentration  be¬ 
tween  lesion  and  surrounding  healthy  tissue,  the  measured 
signal  difference  (or  contrast)  iu  the  scintigram  depends 
strongly  on  imaging  geometry,  since  attenuation  of  gamma  rays 
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exiting  the  lesion,  as.  well  as  the  detector  area  over  which  the 
lesion  image  is  spread,  are  both  functions  of  the  lesipn-to-col- 
limator  distance  [1],  [2],  If  the  location  of  the  lesion  within  the, 
breast  is  known,  it  may  be  possible  to  more  accurately  estimate 
the  relative  radiotracer  concentration  by  making  explicit  cor¬ 
rections  for  attenuation  and  camera  spatial  resolution.  Out  goal, 
is  to  use  the  ability  of  the  X-ray  component  of  the  dual  modality 
system  to  localize  the  lesion  in  the  breast,  thereby  permitting 
these  corrections.  We  present  the  results  of  phantom  studies 
demonstrating  the  feasibility  of  improved  lesion-to-background 
uptake  ratio  assesspaent  via  this  technique.  .  '  ' 

n.  Dual  Modality  System  Description.  .  .  • 

The  digital  X-rhy  and  garoma  ray  detectors  are  mounted  on 
a  con'f/entional  upright  moimnographic  unit  (Lofad  M-ID).  .The 
digital  mammography  detector,  designed  and  developed  by  our 
group  in  collaboration  with  Brandeis  University,  consists,  of  six 
butted  modules  arranged  in  a  2  x  3  may.  Each  module  con¬ 
tains  a  2  k  X  2  k  CCD  bonded  to  the  small  end  of  a  fiber  optic 
taper.  The  large  ends  of  the  six  tapers  are  optically  coupled  to  a 
common  20  cm  X  30  cm  Gd202S:Tb  phosphor.  The  pixel  size  is 
0.046  mm,  and  a  full  digital  mammogram  contains  4224  x.  6264 
pixels,  witii  16  bits  of  intensity  resolution  per  pixel.  Perfor¬ 
mance  of  the  detector  in  terms  of  the  spatial  firequency  depen¬ 
dent  detective  quantum  efficiency  has  been  published  elsewhere 
[3J. 

The  dedicated  gamma  camera  of  the  dual  modality  breast 
imaging  system  consists  of  16  tiled  Hamamatsu  position  sensi¬ 
tive  photomultiplier  tubes  (PSPMTs),  anrariged  in  a  4  x  4'OTay. 
The  PSPMT  array  is  optically  coupled  via  acrylic  light  guides 
to  a  30  X  30  array  ^  Nal(n)  crystals  from  Saint-Gob^.. 
Each  crystal  is  3  mm  x  3  mm  x  6  mm  (thick),  and  the  crystal 
center-to-center  spacing  is  3.3  mm.  Die  energy  resolution  is 
15%  at  140  keV,  Eccaus'e  of  the  compact  design  of  the  PSPMTs 
^d  camera  housing,  the  gainma  camera  active  area  extends 
to  within  8  mm  of  the  patient’s  chest  wall,  a  performance 
comparable  to  that  of  X-ray  mammography  and  far  superior 
to  that  of  clinicd  gamma  cameras  (at  least  several  Gm)^ -Dae 
camera  uses  a  high  resolution  etched  tungsten  collimator 
with  1.25  jnm  square  apertures  (Tecomet  Corp.),  with  height 
L78  cm  (1.73  cm  effective  length  at  140  keV).  The  Collimator 
septa  are  0.2  mm  thick.  The  intrinsic  detector  spktial  resolu¬ 
tion,  averaged  over  the  detector  surface,  is  3.8  mm  FWHM. 
Using  a  — 7%/-|-22%  energy  window,  the  camera  sensitivity  is , 
290  cpm/ yiCi  (efficiency  of  1 .3  x  10“^). 

^'The  current  scaimer,  now  undergoing  clinical  evaluation,  is 
shown  in  Fig.  1.  In  a  typical  dual  modality  study,  the  patient  is 
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Fig.  1.  Dual  modality  scanner  now  tinder  clinical  evaluation  at  the  University 
of  \^ginia.  The  gamma  camera  is  raised  to  the  location  shown  for  the  X-ray 
image,  then  lowered  to  a  position  immediately  above  the  compression  paddle 
for  acquisition  of  the  gamma  image. 

first  positioned  in  a.  manner  similar  to  that  used  in  conventional 
mammography  (i.e.,  with  the  breast  held  between  the  upper  sur-. 
face  of  the  X-ray  detector  and  the  compression  paddle),  except 
with  substantially  reduced  breast  compression.  With  the  gamma  ^ 
camera  pivoted  out  of  the  path  of  the  X-ray  beam  a  digital  mam¬ 
mogram  is  obtained,  and  the  image  appearing  on  the  acquisi¬ 
tion  monitor  is  used  to  verify  that  the  portion  of  the  breast  con¬ 
taining  the  suspicious  lesion  is  approximately  centered  in  the 
field  of  view.  Next,  the  gamma  camera  is  swung  back  into  po¬ 
sition  above  the  breast  and  is  lowered  via  a  linear  translation 
stage  to  a  position  just  above  the  compression  paddle.  Following 
acquisition  of  the  gamma  emission  image  15  minute  ac¬ 
quisition  time),  the  X-ray  and  gamma  ray  images  are  corrected 
for  translational  and  rotational  offsets  and  their  pixel  sizes  are 
rescaled  to  a  common  value  using  a  stored  coregistration  file. 
X-ray  and  gamma  ray  images  are  separately  window/leveled, 
then  the  common  areas  of  each  are  fused  into  a  single  dual 
modality  image.  Details  of  the  procedures  for  merging  the  X-ray 
and  gamma  ray  images  are  essentially  the  same  as  those  used  in 
our  dual  modality  small  animal  scanners  and  are  described  else¬ 
where  [4]. 

In  the  next  generation  of  scanners  now  under  development, 
the  breast  is  supported  by  an  independent  structure  rather  than 
by  the  X-ray  detector.  This  support  structure  is  mounted  on  a 
.  stainless  steel  rod  that  is  attached  to  the  generator  housing  and 
goes  through  the  gantry  arm.  Since  it  is  not  attached  to  the  gantry 
aim,  the  X-ray  and  gamma  ray  detectors  cm  rotate  indepen¬ 
dently  of  the  breast.  Fig.  2  shows  the  gantry  arm  and  support 
structure.  The  gantry  arm  is  driven  by  a  dc  servo  motor  via  a 
reducing  worm  gear,  permitting  rapid  reproducible  angular  dis¬ 
placement  of  the  detectors.  Thus,  it  is  possible  to  obtain  mul¬ 
tiple  views  of  the  breast  in  a  single  configuration  (no  change 


tfrp.  transactions  ON  NITCT.EAR  SCIENCE.  VOL.  50,  NO.  3,  JUNE  2003 


Fig.  2.  Dual  modality  scanner  currently  under  development.  The  breast 
support  is  now  independent  of  the  detectors  and  gantry  arin,  permitting  multiple 
views  of  the  breast  to  b6  ObMietl.  IhB  iJbfVU  lilOLUl  aiid  WuilU  gbai  uatibc-abun 
to  the  right  of  the  vertical  support  frame. 


Fig.  3.  Schematic  diagrmn  of  the  cylindrical  phantom  used  for  the  study.  0^ 
lesion  is  shown,  located  at  a  distance  LCD  from  the  collimator.  The  length  D 
is  the  upper  limit  in  the  integral  of  (1).  The  inner  diameters  (ID)  of  the  lesions 
used  in  the  study  were  1.0  and  0.8  cm.  The  view  shown  is  at  5  =  0° . 

in  compression  orientation)  at  precisely  known  angles.  By  ac¬ 
quiring  two  X-ray  projections  of  the  breast,  the  lesion  can  be 
accurately  localized  within  the  breast,  as  is  done  on  stereotactic 
breast  biopsy  syvStems, 

HI.  Phantom  Studies 

Phantoms  containing  simulated  lesions  were  constructed  to 
measure  the  improvement  in  concentration  ratio  assessment 
through  lesion  localization.  Fig.  3  shows  a  schematic  diagram 
of  one  such  phantom  and  illustrates  the  acquisition  geometry 
used  in  the  study.  Lesions  were  attached  via  low  attenuation 
rods  of  known  length  to  a  central  tube.  Each  rod  screwed 
into  one  of  a  set  of  regularly  spaced  holes  drilled  in  the  tube. 
The  tube  was  coaxial  with  an  external  rod  used  to  mount  the 
phantom  on  a  rotation  stage.  Thus  the  tube’s  axis  coincided 
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Fig.  4.  Image  of  the  phantom  used  to  measure  scatter  fraction  for  the  lesion 
contrast  studies  described. 

with  the  rotation  axis  (the  y  axis  .in  Fig.  3).  The  radial  distance 
from  the  ph^tom  center  (r)  and  the  starting  angular  position 
of  each  lesion  relative  to  the  +z  axis  (cp)  was  precisely  known 
based  on  the  rod  lengths  and  location  of  the  holes  in  the  tube. 
Bends  of  9D°  were  made  in  the  rods  to  minimize  superposition 
of  lesions  and  rods  as  the  phantom  was  rotated. 

For  the  phantom  studies,  a  gamma  camera  very  similar  to  that 
described  above  was  used.  The  P5PMT  ^d  crystal  arrays  were 
nominally  identical,  but  in  this  case  the  collOTator  was  made  of 
lead,  with  hexagonal  holes  with  length  2.1  cm  (2.0  cm  effective 
length  at  140  keV),  and  inner  wall-to-wall  dimension  of  1 .6  mm. 
Using  a  — 8  J%/-|-24%  energ}^  window,  the  gamma  camera  sen¬ 
sitivity  was  230  cpm/jLfcCi  (efficiency  of  1.03  xlO""^),  similar 
to  that  of  the  camera  with  the  tungsten  collimator.  A  phantom 
with  a  cylmdrical  cross  section  was  tested.  The  inner  diameter 
of  the  phantom  was  9.6  cm,  and  its  height  was  10.7  cm..  Lesions 
had  inner  diameters  of  8  mm  and  10  mm,  with  outer  diameters 
of  lb  mm  and  12  mm,  respectively.  Lesions  were  obt^ed  from 
Data  Spectrum  Coip.^  Realistic  lesion-to-background  concen¬ 
tration  ratios  (^  6:1  -10:1)  were  used.  Gamma  images  were 
.  obtained  oyer  a  large  range  of  viewing  angles,  6,  Thus  the  an¬ 
gular  offset  from  the  z  axis  of  the  zdi  lesion  at  wewing  angle  6 
was  (cpi  -f  9),  For  each  view,  and  for  each  lesion,  slices  passing 
flnough  the  lesion  were  taken  from  tjie  image. 

To  quantify  the  effects  of  .attenuation,  point  sources 

were  used  to  directly  measure  the  attenuation  of  the  acrylic  from 
which  the  phantoms  were  constructed.  The  measured  linear  at¬ 
tenuation  was  0.18  cni"^.  The  measured  hnear  attenuation  of 
water  was  0.16  cm"^. 

Scatter  to  primary  ratios  were  determined  empirically  using  a 
ph^tom  the  same  size  as  those  containing  the  simulated  lesions, 
but  containing  a  “cold”  region  (water  only).  An  image  of  the 
ph^tbm  is  shown  in  Fig.  4.  Fig.  5  is  a  plot  of  the  energy  spectra 
obtained  in  the  hot  and  cold  regions  of  the  phantom.  These 
spectra  were  used  to  determine  the  fraction  of  total  counted 
events  arising  from  direct  (not  scattered)  gamma  rays.  For  the 
phantom  data  presented  here,  the'  scatter  fraction  was  0.32. 

Views  were  obtained  by  rotating  the  phantoms  on  a  com¬ 
puter-controlled  rotary  stage  md  holding  the  c^era  stationary. 
Acquisition  times  for  each  view  were  adjusted  according  to  the 
6  hour  half-life  of  so  that  approximately  the  same  total 

counts  were  obtained  for  each  view. 


Scatter  In  Phantom  Measurements 


Channel  Number 


Fig.  5.  Spectra  obtained  using  the  phantom  whose  image  is  shown  in  Fig.  4. 
These  spectra  were ‘used  to  ciculate  the  scatter  fraction  for  a  given  energy 
window.  ^ 

rv.  Theory 

Three  factors  reduce  lesion  coritrast  m  the  images:  1)  attenu¬ 
ation  of  gamma  rays  radiating  from  the  lesion;  2)  partid  volurne 
averaging  (camera  spatial  resolution);  and  3)  scattered  gamma 
rays  whose  detected  energy  falls  writhin  the  energy  window  nsed 
for  acquisition.  We  define  ni  as  the  number  of  epunfr  m  a  lesion 
pixel,  and  as  the  number  of  counts  in  a  background  pD:el. 
Theoretical  expressions  for  ul  ni  that  mclude  the  effects 
of  ganma  camera  spatial  resolution,  attenuation,  and  scatter  are 
obtained'  as  follows.' 

The  total  number  of  counts  in  a  pixel  in  the  image  is  ^veh  by 


71  —  A-w^allP 


(1) 


where  e  is  the  detection  efficiency,  t  is  the  acquisition  time, 
Awaii  is  the  attenuation  due  to  the  acrylic  wall  on  the  camera 
side  of  the  phantom,  p  is  the  pixel  size,  C{z)  is  the  activity  per 
unit  voliune,  D  is  the  thickness  of  the  phantom  above  the  pixel, 
and  iJLyj  is  the  linear  attenuation  coefficient  of  water  for  140  keV . 
gamma  rays.  Note  that  by  using  pyj  only  in  (1),  we  have  ignored 
the  difference  between  the  linear  attenuation  of  water  and  that  of 
the  walls  of  the  lesions,  even  though  some  detected  gainiria  rays 
are  emitted  frpm  background  regions  behind  a  lesion.  This  is 
reasonable  since  the  lesion  walls  are  thin  (^  1  mm)  compared  to 
the  overall  phantom  dimensions  (<^9-15  cm).  Thus,  for  a  pixel 
viewing  the  background  region  of  the  phantom,  relative  to  the 
zero  attenuation  case,  the  reduction  in  counts  due  to  attenuation 
is  .given  by  •  .  ^ 


_  -^wall  [f  ,  exp  (  /i^Z))] 


(2) 


The  total  number  of  image  counts  due  to  gamnia  rays  coming 
from  the  lesion  is  given  by 


total  —  CMeA  (3) 


iffillsborongh,  NC,  USA. 


•where  Ci  is  the  activity  concentration  in  the  lesion,  Vj,  is  the 
lesion  volume,  t  is  the  acquisition  time;  e  is  the  gamma  camera 
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efficiency,  and  A  is  the  attenuation  of  the  water  and  acrylic  lying 
between  the  lesion  and  the  gamma  camera. 

To  take  into  account  the  reduction  in  contrast  due  to  camera 
blur,  we  calculate  the  shape  of  the  lesion’s  projection  onto  the 
detector,  convoluted  with  the  camera’s  two-dimensional  (2-D) 
point  spread  function <PSF).  With  no  blur,  the  2-P  count  density 
functi.on  of  the  image  of  a  spherical  lesion,  projected  onto  a  the 
surface  of  the  detector,  has  the  form 


X  for 


1- 


,./|2 


+  b'r 


1 


■<1 


(4) 


where  R  is  the  lesion  rachus,  xl  and  yi  are  the  x-  and  y-coor- 
dinates  of  the  lesion  center,  x'  =  x-  xl,  t/  =  y  —  Vl,  and  the 
peak  amplitude  JVx,,max  is  set  by  the  condition 


// 


NL(x,y)dxdy  =  total- 


(5) 


In  these  phantom  studies,  the  lesion-to-coUimator  distance 
(LCD)  is  a  function  of  the  phantom  rotation  angle,  as  is  the, 
location,  x,y,  on  the  collimator  surface  onto  which  the  image 
of  the  lesion  projects  (see  Fig.  3).  Including  detector  blur,  lesion 
.  profile  is  given  by  the  convolution  of  the  projected  sphere  profile 
iVjc,  (a:,  r/)  and  the  2-D  Gaussian  PSF 


Nl, y)  =iVi(a:,  y)  ®  g{x,  y) 

.  ,  s  f  ix-XLf  +  {y-yLf\ 

where  g{x,  y)  =  exp  I - j  w 


where  cr  is  equal  to  the  full-width  at  half-maximum  (FWHM) 
of  the  camera  LSF  divided  by  2.355. 

To  find  <T  as  a  fimetion  of  LCD,  the  line  response  function 
(LSF)  of  the  gamma  camera  was  measured  over  a  range  of 
source-to-coUimator  distance  ranging  from  the  minimum  to 
TTiaTiTnnm  LCD  values  used  in  the  phantom  study.  A  thin 
1  rom  inner  diameter)  capillary  tube  filled  with 
solution  was  used  to  produce  a  fine  source  input.  The  capillary 
was  angled  relative  to  the  rows  and  columns  of  the  detector 
crystal  array  so  that  the  phase  relationship  between  the  line 
input  and  the  crystal  sampling  matrix  varied  over  the  length  of 
the  imaged  line;  Camera  resolution  for  each  capillary-to-col- 
limator  distance  was  computed  as  the  FWHM  of  the  image 
of  the  capillary,  averaged  over  its  length.  Fig.  6  shows  the 
measured  FWHM  camera  resolution,  ^ofiles  made  through 
the  line  source  images  indicated  that  the  camera  LSF  is  well 
represented  by  a  Gaussian  function. 

Convolution  with  the  2-D  Gaussian  PSF  broadens  the  shape 
of  y)  with  respect  to  that  of  Nl{x^  y).  The  require¬ 

ment  that 


(7) 
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Fig,  6.  Measured  FWHM  gamma  camera  spatial  resolution  as  a  function  of  the 
source-to-coUimator  distance.  Values  shown  represent  the  resolution  averaged 
over  the  surface  of  the  camera. 

copresponds  physically  to  the  requirement  that  the  total  number 
of  counts  in  the  bluired  profile  equals  that  in  the  unbluired  pro¬ 
file.  This  requirement  sets  the  amplitude  of  NL-;bha{x,  y)- 
The  image  contrast  with  and  without  scatter  radiation  is 

^  tI’Xj  ‘  rift  ^  ^  _  (n.j,  -j-  71-5)  —7  (ri-b  ~i~  n.^) 

nb  ’  *  nb  + ns 

where  n^,  Ub,  refer  to  the  number  of  counts  in  a  pixel  at  the 
center  of  the  lesion  image,  and  in  the  background,  respecfively, 
and  ns  refers  to  the  additional  counts  fiom  scatter.  Here,  we 
make  the  reasonable  assumption  that  the  number  of  scatter 
counts  in  a  lesion  pixel  is  approximately  equal  to  that  in  a 
background  pixel. 

Equation  (8)  indicates  that  the  reduction  in  contrast  due  to 
scattered  radiation  is 

•  Cs 


Ub 


0  rib  Us 


=  i 


nv 


(9) 


nfe  “h 

The  ratio  nsl{nb  -h  n^)  (the  scatter  firaction)  was  determined 
from  the  Spectra  shown  in  Fig.  5,  using  the  energy  window  uti- 
lized  for  image  acquisition.  This  permits  ns  to  be  deterrnined  for 
any  given  value  of  nj  calculated  firom  scatter-firee  conditions. 

The  theoretical  value  for  the  number  of  counts  in  a  back¬ 
ground  pixel  is  obtained  b)'’  using  (1).  The  theoretical  value  for 
the  number  of  counts  in  the  center  pixel  of  the  lesion  image 
is  obtained  by  using  (3)-(7)  to  obtain  the  peak  amplitude  of 
NlMut  {x,y).  This  peak  amplitude  gives  the  quantity  (ni  - nj) . 
Then  (8)  and  (9)  areused  to  obtain  the  calculated  image  contrast. 

V.  Image  Analysis  . 

To  quantify  the  effects  of  attenuation,  spatial  resolution,  and 
scatter  in  the  phantom  image  data,  the  phantom  images  were  an¬ 
alyzed  as  follows.  First,  for  each  lesion,  a  row  of  pixels  going 
through  the  lesion  center  was  extracted.  Next  several  rows  near 
the  first  row,  but  containing  only  background  were  extracted. 
The  background-only  rows  were  averaged  to  obtain  a  back¬ 
ground-only  row  in  the  region  of  the  lesion.  This  averaged  back- 
ground-only  row  was  then  subtracted  firom  the  row  containing 
tile  lesion  on  a  pixel-by-pixel  basis.  The  resulting  difference 
row  contains  only  the  profile  of  the  lesion  and  noise.  For  sim¬ 
plicity,  a  Gaussian  function  was  fit  to  the  profile  to  obtain  a 
noise-free  estimate  of  the  peak  amplitude.  For  the  (small)  lesion 
sizes  and  (relatively  large)  lesion-to-collimator  distances  used  in 
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TABLE  I 

Parameters  for  Phantom  Image  Acquisition 


Quantity 

Symbol 

Value 

Camera 

detection 

efficiency 

e 

diinensionless 

1.03x10^ 

Cb 

jiCi/cm^ 

1.4 

Lesion  activity 
concentration 

Cl 

pCi/cm^ 

14.5 

Acquisition 
time  per  view 

t 

1200  seconds  for 
first  vievv,  decay 
conected 
thereafter 

Pixel  area 

Ap  • 

cm^ 

0.11 

Scatter 

fraction 

SF 

diTTi.ensionle.ss 

0.32  • 

LCD 

cm 

-  4  (minimum) 

-  9  (maximum) 

Lesion  Profile  Fitting 


Fig.  7.  Plots  of  experimentally  determined  lesion  profile  (straight  line 
segments),  and  the  function  Pi,  (smooth  curve),  given  by  the  convolution  of 
(3).  ■  ■ 

this  study,  a  Gaussian  function  was  a  reasonable  fit  to  the  rdea- 
sured  lesion  profile.  The  peak  amplitude  of  the  fit,  along  with 
the  pixel  value  in  the  corresponding  pixel  of  the  averaged  back¬ 
ground-only  row,  were  used  to  calculate  the  measured  image 
contrast.  ■ '  * 

VI.  Results 

Table  I  summarizes  the  parameter  values  used  in  the  acqui¬ 
sition.  Fig.  7  shows  an  example  of  a  profile  through  a  lesion 
phantoin,  along  with  a  slice  through  the  center  of  the  convolu¬ 
tion  defined  in  (6).  The  particular  example  shown  is  for  R  = 
0.5  cm  and  cr  =  0.27  cm.  Although  the  relatively  small  lesion 
size  causes  its  image  to  be  spread  over  a  small  number  of  pixels, 
the  figure  shows  that  the  lesion  profile  is  well  fit  by  the  expres¬ 
sion  of  (6).  Figs.  8-10  are  plots  of  the  measured  and  calculated 
values  for  nx,  +  tTs,  and  the  contrast,  for  a  study  using 

the  cylindrical  phantom.  Views  were  obtained  at  15°  intervals 
over  345°. 

The  errors  in  the  predicted  values,  relative  to  the  measured 
values  are  given  in  Table  n. 


Fig.  8,  Measured  and  calculated  values  of  near  one  of  the  lesions 

in  the  cylindrical  phantom.  Views  were  obtained  over  345^,  with  15®  between 
views.  The  minimum  LCD,  at  0  =  0® ,  was  3.85  cm,  and  the  maximum  LCD, 
■at  ^  =  180®,  was  9.15  cm. 


Fig.  9.  Measured  and  calculated  values  of  ux,  -H  for  one  of  the  lesions 
in  the  cylindrical  phantom.  Views  were  obtained  over  345®,  with  15®  between 
views.  The  minimum  LCD,  at  0  =  0® ,  was  3.85  cm,  and  the  maximum  LCD, 
at  ^  =  180®,  was  9.15  era.  ■ 


Fig.  10.  Measured  and  calculated  values  of  lesion  image  contrast  for  the  lesion 
whose  data  are  shown  in  Figs.  8  and  9.  Views  were  obtained  over  345®  s,  with 
15®  between  views.' 

vn.  Discussion  and  Conclusion 

The  above  results'  indicate  that  the  theoretical  construct 
models  the  true  effects  from  attenuation,  spatial  resolution,  and 
scatter  with  reasonable  accuracy.  In  these  studies,  the  lesion 
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TABLE  n 

AVERAGE  Error  and  RMS  Deviation  in  the  Error,  Over  345° ,  of  • 
Tffi  Quantthes  Plotted  in  Figs.  8~10 


location,  lesion  size,  attenuating  thickness  of  overlying  breast 
tissue,  lesion-to-coUimatof  distance,  and  scatter-to-piimary 
ratio  were  known  d  priori  from  known  phantom  dimenMons 
and  direct  measurement  of  scatter.  In  a  climcal  context,  tiie 
breast  shape  will  be  determined  from  the  X-ray  images  and 
the  compressed  breast  thickness.  The  lesion  size  can  be  esti¬ 
mated  based  on  both  the  X-ray  and  gamma  ray  images.  The 
scatter-to-primary  ratio  is  a  function  of  the  location  in  the  breast 
and  increases  with  decreasing  distance  to  the  chest  wall.  .In 
our  current  dual  modality  breast  scanner,  we  routinely  compile 
energy  spectra  for  several  separate  contiguous  regions  going 
between  the  chest  wall  and  anterior  edges  of  the  detedtor.  These 
spectra  could'be  used  to  estimate  the  scatter-to-primary  ratio  at 
a  given  location  within  the  breast.  Thus,  if  the  location  within 
the  breast  of  the  lesion  is  known,  these  formalisms  can  be  used 
to  correct  the  measured  (image)  contrast  so  that  it  is  closer 
to  the  true  (radioisotope)  contrast.  The  use  of  multiple  view 
X-ray  imaging  in  the  dual  modality  system  vdll  permit  this 
localization  to  be  made  prior  to  gamma  image  acquisition.  It  is 
anticipated  that  such  techniques  may  improve  quantification  of 
radiotracer  uptake  m.  scintimammography. 
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There  are  several  challenges  to  clinical  implementation  of 
this  approach.  First,  there  will  be  a  certain  ntoiber  of  cases  in 
which  localization  of  the  lesion  and  estimation  of  its  size  via 
X-ray  may  be  dhBficulthecause  of  surrounding  radiodense  breast 
tissue.  One  possible  solution  for  such  cases  is  to  obtain  lesion 
position  information  from  stereoscopic  gamma  images.  Second, 
regional  variations  in  scatter-to-primary  ratio  within  the  breast 
may  make  estimation  of  the  scatter  fraction  in  the  region  of  the 
lesion  difficult.  We  are  currently  assessing  the  magnitude  and 
nature  of  these  variations  in  data  obtained  in  our  ongoing  planar 
du^  modality  clinicai  study.  Third,  the  background  radiotracer 
is  likely  to  vary  within  the  breast,  due  to  normal  differences  in 
uptake  between  various  tissues.  This  will  make  assessment  of 
the  contrast  dependent'on  the  region  of  the  image  chosen  for 
measurement  of  the  quantity  -f  . 


References 

[1]  I.  N.  Weinberg,  R.  Pani,  R-  PcUegrini,  K  Scopinaro,  G.  De  Vincentis,  A. 
Pergola,  and  A.  Soluri,  “Small  lesion  visibility  in  sointimammography,” 
IEEE  Trans.  Nucl  Set,  vol.  44,  pp.  1398-1402, 1997. 

[2]  S.  Majewski,  D.  Kieper,  E.  Curran,  C.  Keppel,  B.  Kross,  A.  Palumbo;  V. 
Popov,  A  G.  Weisenberger,  R.  Welch,  R.  Wqjcik,'M.  B.  WiUiams,  M. 
More,  and  G.  Zhang,  “Optimization  of  dedicated  scintimainmogra|)hy 
procedure  using  detector  prototypes  and  eompr^sible  phantoms,”  ZEES 
Trans.  Nucl.  Set,  vol.  3,  pp.  822—829,  2001. 

[3]  M,  B.  Williams,  P.  U.  Simoni,  L.  Smilowitz,  M.  Stanton,  and  W.  Phillips, 
“Analysis  of  the  detective  quantum  efficiency  of  a  developmental  de¬ 
tector  for  digital  Tnammography,”  Med.  Phys.,  vol.  26,  pp.  2273—2285, 
1999. 

[4]  M.  B.  Williams,  V.  Gaibis-Reig,  A.  Goode,  P.  Simoni,  S.  Majewski,  A. 

Weisenberger,  R.  Wojcik,  W.  Phillips,  and  M.  Stanton,  “Dual  modality 
system  for  small  animal  imagirig,”  RSNA  Electron.  vol.  3, 1999. 


Paper  #M3-241 


appendix  2 


Analysis  of  Position-Dependent  Compton  Scatter  in 
Scintimammography  With  Mild  Compression 


Mark  B.  Williams,  Member,  IEEE,  Deepa  Narayanan,  Mitali  J.  More,  Patricia  J.  Goodale,  Stan  Majewski, 

Douglas  Kieper 


Abstract—  In  breast  scintigraphy  using  ^^"^c-sestamibi  the 
relatively  low  radiotracer  uptake  in  the  breast  compared  to  that 
in  other  organs  such  as  the  heart  results  in  a  large  fraction  of  the 
detected  events  being  Compton  scattered  gamma  rays.  In  this 
study,  our  goal  was  to  determine  whether  generalized  conclusions 
regarding  scatter-to-primary  ratios  at  various  locations  within 
the  breast  image  are  possible,  and  if  so,  to  use  them  to  make 
explicit  scatter  corrections  to  the  breast  scintigrams.  Energy 
spectra  were  obtained  from  patient  scans  for  contiguous  regions 
of  interest  (ROIs)  centered  left  to  right  within  the  image  of  the 
breast,  and  extending  from  the  chest  wall  edge  of  the  image  to  the 
anterior  edge.  An  anthropomorphic  torso  phantom  with  fillable 
internal  organs  and  a  compressed-shape  breast  containing  water 
only  was  used  to  obtain  realistic  position-dependent  scatter-only 
spectra.  For  each  ROI,  the  measured  patient  energy  spectrum 
was  fitted  with  a  linear  combination  of  the  scatter-only  spectrum 
from  the  anthropomorphic  phantom  and  the  scatter-free 
spectrum  from  a  point  source.  We  found  that  although  there  is  a 
very  strong  dependence  on  location  within  the  breast  of  the 
scatter-to-primary  ratio,  the  spectra  are  well  modeled  by  a  linear 
combination  of  position-dependent  scatter-only  spectra  and  a 
position-independent  scatter-free  spectrum,  resulting  in  a  set  of 
position-dependent  correction  factors.  These  correction  factors 
can  be  used  along  with  measured  emission  spectra  from  a  given 
breast  to  correct  for  the  Compton  scatter  in  the  scintigrams. 
However,  the  large  variation  among  patients  in  the  magnitude  of 
the  position-dependent  scatter  makes  the  success  of  universal 
correction  approaches  unlikely.^ 

Index  Terms —  breast  imaging,  Compton  scatter, 

multimodality,  scintimammography 

1.  Introduction 

During  the  past  several  years,  we  have  been  developing  a 
breast  scanning  system  that  combines  digital  x-ray 
mammography  with  breast  scintigraphy  using  a  dedicated 
small  field  of  view  gamma  camera  [1].  The  x-ray  and  gamma 
ray  detectors  are  mounted  on '  an  upright  mammography 
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gantry,  and  the  x-ray  transmission  and  gamma  emission  views 
are  obtained  sequentially  while  the  breast  is  held  in  a  fixed 
configuration  using  mild  compression.  The  chest  wall  edges 
of  the  x-ray  and  gamma  ray  detectors  coincide,  and  the 
gamma  camera  is  centered  left-to-right  within  the  field  of 
view  of  the  x-ray  detector.  The  fields  of  view  of  the  gamma 
camera  and  x-ray  detector  are  10  cm  x  10  cm  and  19  cm  x  28 
cm,  respectively.  In  an  ongoing  clinical  evaluation,  the 
resulting  dual  modality  images  have  proven  effective  for 
differentiation  of  benign  and  malignant  breast  masses, 
compared  to  prone  scintimammography  alone,  with  our 
dedicated  camera. 

The  presence  of  counts  in  the  gamma  image  due  to 
Compton  scattered  gamma  rays  can  reduce  lesion  contrast  [2]. 
In  this  study,  our  goal  was  to  determine  whether  generalized 
conclusions  regarding  scatter-to-primary  ratios  at  various 
locations  within  the  image  are  possible,  and  if  so,  to  use  them 
to  make  explicit  scatter  corrections.  In  a  typical  gamma  image 
acquisition  using  the  dual  modality  system,  a  two-dimensional 
image  histogram  is  updated  each  time  a  gamma  interaction 
occurs  that  deposits  an  amount  of  energy  in  the  detector  that 
falls  w  ithin  a  predetermined  e  nergy  w  indow.  H  owever,  i  n  a 
number  of  recent  patient  studies,  gamma  emission  data  were 
acquired  in  list  mode  rather  than  in  the  conventional  single 
histogram  mode.  In  list  mode,  the  location  on  the  detector 
surface  of  the  g  amma  r  ay  i  nteraction,  t  he  t  ime  a  t  w  hich  t  he 
interaction  occurred,  and  the  total  energy  of  the  interaction  are 
stored  individually  for  each  detected  gamma  event.  The 
resulting  list  of  x,  y,  t,  and  E  values  can  then  be  replayed  with 
any  desired  energy  discrimination  criteria.  In  addition,  energy 
spectra  corresponding  to  any  arbitrary  region  of  the  detector 
surface  may  be  obtained.  We  have  used  this  capability  to 
study  the  spatial  dependence  of  the  Compton  scatter  coming 
from  the  breast. 

II.  Methods 

A.  Imaging  System 

The  dual  modality  breast  imaging  system  at  the  University  of 
Virginia  consists  of  a  digital  x-ray  detector  and  a  dedicated 
gamma  camera  mounted  on  a  single  upright  unit  (Lorad  Mill). 
The  x-ray  detector,  built  at  Brandeis  University,  Waltham, 
MA,  consists  of  six  butted  modules  arranged  in  a  2  x  3  array. 
Each  module  contains  a  2k  x  2k  CCD  bonded  to  the  small  end 
of  a  3.3:1  fiber  optic  taper.  The  large  ends  of  the  six  tapers  are 
optically  coupled  to  a  common  20  cm  x  30  cm  Gd202S:Tb 
phosphor.  The  pixel  size  is  46  microns,  and  a  full  digital 
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mammogram  contains  4224  x  6264  pixels,  with  16  bits  per 
pixel.  The  inactive  border  at  the  chest  wall  edge  is  7  mm  wide. 

The  dedicated  gamma  camera  of  the  dual  modality  breast 
imaging  system,  built  at  the  Thomas  Jefferson  National 
Accelerator  Facility  (TJNAF),  is  mounted  on  the  same  gantry 
and  is  moved  out  of  the  way  while  acquiring  the  x-ray  image. 
The  gamma  camera  consists  of  a  4  x  4  array  of  one  inch 
square  position  sensitive  photo  multiplier  tubes  (PSPMTs) 
(Hamamatsu  R7600-00-C8).  The  PSPMT  array  is  optically 
coupled  via  acrylic  light  guides  to  a  30  x  30  array  of  NaI(Tl) 
crystals.  Each  crystal  is  3.0  mm  x  3.0  mm  x  6  mm  thick,  and 
the  crystal  center-to-center  spacing  is  3.3  mm.  Because  of  the 
compact  design  of  the  PSPMTs  and  camera  housing,  the 
gamma  camera  is  able  to  image  to  within  8  mm  of  the 
patient’s  chest  wall.  The  camera  uses  a  high-resolution  etched 
tungsten  parallel  hole  collimator  with  square  apertures.  The 
single  crystal  energy  resolution  is  ~15%  at  140  keV.  Images 
are  typically  obtained  using  a  +22/-7%  energy  window  set 
around  the  140  keV  photopeak  of  ^^”"Tc. 

B.  Patient  Imaging 

During  the  past  three  years,  the  dual  modality  system  has 
been  used  in  a  clinical  study  conducted  at  the  Wmary  Care 
Center  at  the  University  of  Virginia  to  test  its  ability  to 
distinguish  between  malignant  and  benign  lesions  [1]. 
Participants  are  recruited  from  patients  who  have  already 
undergone  screening  mammograms  and  are  scheduled  for 
biopsy. 

For  this  study,  the  patient  is  given  an  intravenous 
injection  of  25-30  mCi  of  ^^'"Tc-sestamibi  in  the  arm 
contralateral  to  the  suspicious  breast.  The  suspicious  breast  is 
then  positioned  on  the  x-ray  detector.  The  compression  paddle 
provides  mild  compression  to  maintain  a  consistent  breast 
configuration  for  both  imaging  modalities.  An  x-ray  image  is 
obtained  to  confirm  that  the  suspicious  lesion  will  be  within 


Fig.  1.  Example  patient  gamma  image  showing  five  ROIs,  going  from 
the  chest  wall  edge  of  the  image  to  the  anterior  edge,  along  the  midline 
of  the  breast. 
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Fig.  2d.  Energy  spectrum  from  ROI 4  of  Fig.  1. 
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Fig,  2e.  Energy  spectrum  from  ROI  5  of  Fig.  1.  There  are  almost  no 
primary  counts  because  this  ROI  is  anterior  to  the  breast. 

the  10  cm  X  10  cm  field  of  view  of  the  gamma  detector  when 
it  is  swung  into  position.  Gamma  image  acquisition  is  begun 
once  correct  location  of  the  lesion  in  the  x-ray  image  is 
verified.  The  typical  acquisition  time  is  --1 0-1 5  minutes.  The 
above  procedure  (x-ray  and  gamma  ray  imaging)  is  then 
performed  on  the  non-suspicious  breast.  For  the  studies 
described  here,  acquisition  was  done  in  list  mode,  whereby  the 
time,  position,  and  energy  of  each  gamma  interaction  in  the 
detector  is  stored.  The  threshold  setting  for  the  discriminator 
gating  the  ADCs  was  set  just  above  the  electronic  noise,  so 
that  the  entire  energy  spectrum  was  recorded. 

C.  Image  Analysis 

Individual  energy  spectra  were  obtained  for  each  of  five 
contiguous  3.3  cmx  2  cm  regions  of  interest  (ROIs)  in  the 
gamma  images  of  both  breasts  of  each  patient.  Die  ROIs  were 
centered  left  to  right  in  the  image,  and  extended  from  the  chest 
wall  edge  of  the  image  to  the  anterior  edge  (see  Fig.  1).  ROI  5 
is  the  region  farthest  from  the  chest  wall  edge  while  ROI  1  is 
the  region  closest  to  the  chest  wall  edge.  Fig.  2(a)-(e)  shows 
the  energy  spectra  for  each  of  the  five  ROIs  are  shown  in  Figs. 


Fig.  3.  Photograph  of  the  anthropomorphic  torso  phantom,  being  imaged 
by  the  dual  modality  scanner.  The  gamma  camera,  above  the  5  cm 
simulated  compressed  breast,  is  obtaining  a  scatter-only  image  because 
the  breast  was  filled  with  water  but  no  ®®^Tc. 

2a)-e).  The  horizontal  axes  of  the  spectra  are  normalized  so 
that  channel  200  corresponds  to  140  keV.  As  we  go  from  ROI 
5  to  1,  in  the  direction  towards  the  chest  Wall,  we  see  that  the 
height  of  scatter  peak  increases  relative  to  that  of  the  primary 
peak.  The  spectrum  from  ROI  5  has  few  primary  (photopeak) 
counts  because  in  this  example  the  breast  does  not  extend  into 
that  part  of  the  image.  It  is  apparent  from  the  spectra  that  even 
with  a  -7%  1  ower  e  nergy  1  evel,  some  s  catter  is  1  ikely  to  be 
present  in  the  image,  particularly  in  regions  near  the  chest 
wall.  Overlap  between  photopeak  counts  and  scatter  counts 
makes  it  very  difficult  to  directly  correct  for  the  unwanted 
scatter  counts.  Therefore,  our  objective  was  to  decompose  the 
spectra  into  their  primary  and  scatter  components,  so  that  the 
scatter  fraction  could  be  estimated. 

D.  Emission  Spectra  Decomposition 

We  model  the  clinical  breast  emission  spectra  as  a  linear 
combination  of  a  scatter-free  spectrum  and  a  scatter-only 
spectrum.  That  is,  for  the  i**'  ROI,  the  emission  spectrum  is: 

=  ^0  /  Sp  +  a^  ^  ^  (1) 

where  ao.i  and  aij  are  spatially  varying  coefficients  that  scale 
the  scatter-fi*ee  spectrum,  Sp,  and  scatter-only  spectrum,  Ss,i 
respectively. 

The  scatter-ffee  spectrum  is  that  from  a  point  source  of 
^^'"Tc,  and  is  independent  of  position  in  the  image. 

The  scatter-only  spectrum  varies  strongly  from  one  location 
to  another,  as  can  be  readily  seen  in  Figs.  2a)-e).  To  obtain 
estimates  of  the  scatter-only  spectra,  we  used  an 
anthropomorphic  torso  phantom  with  individually  fillable 
internal  organs  and  torso  (Radiological  Support  Devices, 
Long  Beach,  CA).  For  this  study,  externally  attached  breasts 
with  uniform  thickness  (5  cm)  in  the  cranio-caudal  direction 
were  used  to  s imulate  c  ompressed  breasts.  The  radioactivity 
concentrations  used  in  the  various  phantom  regions  are  given 
in  Table  I.  The  activity  concentration  distribution  used  was 
based  on  the  bio-distribution  of  ^^"’Tc-sestamibi  at  5  minutes 
post  injection  [3]. 


Fig.  4a.  Energy  spectrum  from  ROI  1,  (Patient  spectrum)  along  with  the 
scaled  primary  and  scatter-only  spectra,  and  their  sum  (Fitted  spectrum). 
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Fig.  4b.  Energy  spectrum  from  ROI  2,  along  with  the  scaled  primary  and 
scatter-only  spectra,  and  their  sum. 
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Fig.  4c.  Energy  spectrum  from  ROI  3,  along  with  the  scaled  primary  and 
scatter-only  spectra,  and  their  sum. 


The  anthropomorphic  torso  phantom  was  scanned  using  the 
dual  modality  breast  scanner  and  the  same  positioning  as  was 
used  for  the  clinical  scans,  although  no  x-ray  images  were 
obtained.  Fig.  3  is  a  photograph  of  the  setup.  Emission  spectra 
were  obtained  from  the  anthropomorphic  phantom  breasts 
using  the  same  ROIs  as  shown  in  Fig.  1.  Using  these  spectra 
as  the  Ss,i  of  (1)  and  the  spectra  obtained  from  the  point  source 
as  Sp,  the  coefficients  aQ  and  aj  were  found  for  each  of  the 
five  ROIs  using  a  least-squares  curve  fitting  routine  that 
minimizes  the  error  between  the  expression  in  (1)  and  the 
measured  clinical  emission  spectrum,  Sg.  The  criterion  used  to 


establish  goodness  of  fit  was  that  the  relative  decrease  in  chi- 
squared  value  from  one  iteration  of  the  curve  fitting  routine  to 
the  next  was  less  than  1.1x10"^.  Figs.  4a)  -  4e)  show  the  fits  to 
the  example  clinical  spectra  of  Fig.  2.  For  each  of  the  five 
ROIs,  the  patient  spectrum,  the  fitted  spectrum  S,  and  its 
components  aoSp  and  aiSj  are  shown. 

The  position-dependent  scatter  fraction,  defined  as  the 
ratio  of  scatter  counts  to  total  counts,  was  then  obtained  as 
follows.  For  a  selected  energy  window,  the  number  of  counts 
in  the  scaled  scatter-only  spectrum  (aiSg)  falling  within  the 
energy  window  was  determined.  The  ratio  of  this  number  to 
the  total  number  of  counts  in  the  original  patient  spectrum 
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falling  within  the  energy  window  is  the  scatter  fraction.  Fig. 
4f  shows  an  expanded  view  of  the  +22/-7%  energy  window, 
with  the  portions  of  aiSj  and  Sc  that  fall  within  it  in  ROI  4. 
The  scatter  fraction  is  the  ratio  of  the  area  under  the  aiSj  curve 
to  that  under  the  Sc  curve. 

Once  the  scatter  fractions  for  each  ROI  have  been 
calculated,  the  images  are  scatter  corrected  as  follows.  The 
scatter  fraction  for  each  ROI  i  s  attributed  to  the  pixel  at  its 
center.  Linear  interpolation  in  the  chest-anterior  direction  is 
used  to  estimate  scatter  fractions  for  the  other  pixels.  We 
assume  a  constant  scatter  fraction  in  the  left-right  direction 
(see  Fig.  7).  Each  pixel  is  scatter  corrected  by  subtracting  a 


number  of  counts  equal  to  its  uncorrected  count  times  the 
scatter  fraction  for  that  pixel. 

III.  Results 

We  have  characterized  the  spatial  dependence  of  the 
scatter  fraction  in  eleven  cases  (twenty-two  dual  modality 
breast  scans).  Figs.  5  and  6  are  plots  of  the  scatter  fraction 
versus  the  distance  in  centimeters  from  the  chest  wall,  for 
cranio-caudal  views  of  the  left  and  right  breasts,  respectively, 
of  seven  patients.  The  ROIs  in  which  the  scatter  fractions 
were  measured  were  2.0  cm  x  3.3  cm  (6  x  10  pixels),  and  were 
placed  along  the  midline  of  the  breast.  The  plots  show  that  for 
each  case,  there  is  a  distinct  decrease  in  the  scatter  fraction 
with  increasing  distance  from  the  chest  wall.  In  the  case  of 
patient  52,  the  scatter  fraction  increases  sharply  beyond  6.5 
cm  from  the  chest  wall  because  in  that  particular  case  the 
compressed  breast  only  extended  -6  cm  from  the  chest  wall. 
The  figures  also  demonstrate  that,  although  the  scatter  fraction 


Distance  from  Chest  Wail  (cm) 


Distance  from  Chest  Wall  (cm) 

Figs.  5  and  6.  Scatter  fraction  versus  distance  from  the  chest  wall  in  the 
left  and  the  right  breast.  Curves  are  labeled  according  to  case  number,  with 
data  points  connected  by  lines  for  clarity. 
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Fig.  7:  Composite  of  5  spectra  taken  from  5  ROls  arranged  along  a  line  in 
the  medial-lateral  direction.  Channel  100  corresponds  to  140  keV. 

decreases  with  increasing  distance  from  the  chest  wall  for  all 
patients,  there  is  a  large  variation  among  patients  in  the  value 
of  the  scatter  fraction  at  a  given  distance. 

We  also  examined  the  change  in  scatter  fraction  within 
the  breast  going  in  the  left-right  (or  medial-lateral)  direction. 
We  found  that  the  scatter  fraction  was  essentially  constant. 
Fig,  7  illustrates  the  consistency  of  the  shape  of  the  emission 
energy  spectrum  along  this  direction  in  an  example  image. 
The  figure  shows  a  composite  of  5  spectra  taken  from  5  ROIs 
arranged  along  a  line  in  the  medial-lateral  direction.  They  are 
virtually  indistinguishable.  We  found  that  this  was  true  for 
each  of  the  breast  images  we  tested,  unless  there  was  an 
abnormality  present,  such  as  a  lesion  with  high  radiotracer 
uptake. 

We  evaluated  the  dependence  of  the  Compton  scatter  on 
breast  thickness.  Figs.  8  and  9  show  the  scatter  fraction  versus 
compressed  breast  thickness,  for  the  left  and  right  breasts, 
respectively.  The  scatter  fractions  shown  were  measured  in  a  1 
cm  X  1  cm  (3  X  3  pixel)  ROI  centered  3.5  cm  from  the 
patient’s  chest  wall.  Each  plot  shows  a  trend  towards 
increasing  scatter  with  increasing  breast  thickness,  however 
there  is  one  case  (compressed  thickness  of  5.5  cm  for  both 
breasts)  in  which  exceptionally  high  scatter  fractions  were 
obtained.  The  reason  for  this  outlier  is  unknown  at  present. 

Fig.  10  shows  an  example  of  the  result  of  scatter 
correction.  The  figure  contains  two  profiles,  e ach  c onsisting 
of  the  same  row  of  pixels,  running  in  the  chest  to  anterior 
direction.  One  profile  is  the  uncorrected  profile  from  the 
original  scintigram,  and  the  other  is  the  same  profile,  but  after 
removal  of  scatter  counts.  As  can  be  predicted  from  Figs.  5 
and  6,  the  number  of  scatter  counts  subtracted  increases  with 
decreasing  distance  from  the  chest  wall.  The  region  of 
uniform  breast  thickness  (where  the  breast  is  in  contact  with 
the  compression  paddle)  corresponds  to  the  range  from 
approximately  2  -  8  cm  from  the  chest  wall.  The  pixel  count 
in  this  region  is  noticeably  more  constant  following  scatter 


Fig.  8:  Scatter  fraction  versus  compressed  breast  thickness  3.5  cm  from 
the  chest  wall,  left  breast. 


correction,  as  would  be  expected  if  the  counts  were  mainly 


primary  counts  arising  from  a  volume  of  uniform  thickness 


and  concentration. 


IV.  Discussion  AND  Conclusions 

Although  we  have  a  limited  number  of  cases  thus  far, 
several  conclusions  can  already  be  made.  In  all  cases  studied, 
there  is  a  marked  decrease  in  Compton  scatter  with  increasing 
distance  from  the  chest  wall.  However,  the  value  of  the  scatter 
fraction  at  a  given  distance  from  the  chest  wall  varies 
substantially  between  patients  (see  Figs.  5  and  6).  This 
variation  indicates  that  ‘universal’  scatter  corrections  based 
simply  on  distance  from  the  chest  wall  are  not  likely  to  work. 
Even  after  correction  for  variation  in  breast  size,  by  examining 
scatter  fraction  at  a  given  fraction  of  the  chest  wall-nipple 
distance,  we  found  a  similarly  large  inter-patient  variation. 

Although  the  amount  of  scatter  radiation  varies  strongly  in 
the  chest-anterior  direction,  we  found  very  little  variation  in 
the  medial-lateral  direction.  In  addition,  the  scatter  fractions  at 
a  given  location  in  the  left  and  right  breasts  of  any  one  patient 
are  quite  similar.  This  invariance  was  somewhat  unexpected, 
given  the  strong  uptake  of  sestamibi  by  the  myocardium. 

We  also  examined  our  dual  modality  images  to  determine 
whether  there  was  any  correlation  between  radiographic 
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Fig.  10.  Profile  (single  row  of  pixels)  running  in  the  chest  wall  to  anterior 
direction,  and  located  at  the  breast  midline.  Profiles  are  shown  before  and 
after  scatter  counts  have  been  removed.  The  region  of  the  image 
corresponding  to  the  portion  of  the  breast  with  uniform  thickness  is  seen  to 
have  more  uniform  pixel  count  following  scatter  correction. 

(mammographic)  density  and  the  scatter  fraction.  No  apparent 
dependence  on  radiographic  density  was  observed. 

In  general,  the  scatter  fraction  at  a  given  distance  from  the 
chest  wall  increases  with  increasing  breast  thickness. 
However,  the  results  shown  in  Figs.  8  and  9  do  not  take  into 
account  variation  among  patients  in  radiotracer  concentration 
in  the  various  parts  of  the  body.  The  ratio  of  breast  uptake  to 
myocardial  uptake,  for  example,  could  have  a  significant 
impact  on  the  measured  scatter  fraction.  While  we  have  no 
way  of  measuring  uptake  in  the  heart  and  liver,  an  estimate 
can  be  made  of  the  average  radioactivity  concentration  in  the 
breast  in  the  region  in  contact  with  the  compression  paddle. 
Using  the  known  gamma  camera  detection  efficiency  (^lO"^), 
the  measured  compressed  thickness,  and  assuming  an 
attenuation  for  breast  tissue  approximately  equal  to  that  of 
water  (p  0.155  cm"^),  we  calculated  the  mean  breast 
radioactivity  concentration  at  the  center  of  ROI  2.  The  mean 
value  obtained  from  12  images,  6  patient  studies,  is  0.2041 
pCi/cm^,  with  an  inter-patient  standard  deviation  of  0.0672 
pCi/cm^.  This  fairly  large  (33%  of  the  mean)  RMS  variation 
among  patients  in  concentration  could  partially  explain  the 
large  variation  we  observe  in  scatter  fraction  at  any  given 
distance  from  the  chest  wall. 

Position-dependent  removal  of  scatter  counts  via  estimation  of 
the  scatter  fraction  has  the  anticipated  effect  of  making  the 
pixel  count  in  the  uniform  thickness  region  of  a  normal  breast 
more  uniform.  However,  even  after  removal  of  scatter  counts, 
we  observe  that  the  total  pixel  count  still  increases  slightly  in 
the  region  near  the  chest  wall.  This  could  be  the  result  of 
incomplete  scatter  removal,  or  may  simply  be  indicative  of  a 
typical  pattern  of  variation  in  radiotracer  concentration 


because  of  normal  breast  physiology  (e.g.  greater  vascular 
density  near  the  chest  wall,  or  the  proximity  of  muscle). 
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ORGAN 

VOLUME 

(ml) 

Heart  (walls) 

295 

1.0 

Liver 

1015 

3.9 

Lung  (right) 

1162 

0.51 

Lung(left) 

882 

0.51 

Torso 

6500 

0.33 

Breast  (right) 

566 

0 

Breast  (left) 

566 

0 

Table  1:  Concentrations  of  radioactivity  used  in  the  anthropomorphic 
phantom  to  obtain  scatter-only  breast  emission  spectra. 
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X-ray  Stereotactic  Lesion  Localization  in 
Conjunction  with  Dedicated  Scintimammography 
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B.  Williams,  Member,  IEEE 


Abstract--  We  are  developing  a  dual  modality  system  that 
combines  digital  x-ray  mammography  with  gamma  emission 
scintigraphy  on  an  upright  mammography  gantry.  The  breast  is 
held  under  mild  compression  by  a  support  structure  that  is 
independent  of  the  detectors.  The  x-ray  source  and  detectors  can 
be  rotated  around  a  fixed  rotation  axis  permitting  multiple  views 
of  the  breast  with  fixed  compression.  Two  such  views  can  be 
combined  as  a  stereotactic  pair  to  obtain  the  three-dimensional 
location  of  breast  lesions.  Information  about  the  location  of  the 
lesion  within  the  breast  permits  corrections  for  attenuation  and 
detector  spatial  resolution,  resulting  in  more  accurate  estimation 
of  the  true  lesion-to-background  concentration  ratio,  based  on 
the  image  lesion-to-background  counts  ratio.  In  this  paper,  we 
describe  the  model  used  to  make  these  corrections,  and  present 
the  results  of  the  phantom  experiments  designed  to  test  the 
accuracy  of  our  calculations. 

Index  Terms —  breast  imaging,  technetium-99m  sestamibi 
uptake,  multimodality,  scintimammography,  stereotactic 
localization 

I.  Introduction 

Accurate  preoperative  staging  of  breast  cancer  is  required 
to  plan  surgery  and  postoperative  therapy.  Scopinaro  et 
al.  have  suggested  that  the  uptake  of  technetium-99m 
sestamibi  is  correlated  with  angiogenesis  and  thus  can  be  a 
potential  marker  for  breast  cancer  invasiveness  [1].  Other 
researchers  have  examined  the  degree  of  correlation  between 
sestamibi  uptake  and  the  traditional  indicators  of  malignancy 
such  as  desmoplasia  and  mitotic  activity  [2].  All  these  studies 
require  quantification  of  the  differential  radiotracer  uptake 
between  the  lesion  and  the  surrounding  healthy  breast  tissue. 
Since  the  depth  of  the  lesion  plays  an  important  role  in 
determining  the  contrast  of  the  lesion  in  the  image  [3],  it  must 
be  taken  into  account  if  accurate  uptake  values  are  to  be 
obtained.  The  attenuation  due  to  the  structure  overlying  the 
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Fig.  1.  Image  showing  2  lesions  of  equal  size  and  equal  radioactivity 
concentration  but  placed  at  varying  depth  from  the  detector  surface.  The 
image  contrast  for  lesion  1  is  0.26  and  lesion  2  is  0.42.  The  true 
concentration  ratio  between  the  lesions  and  the  background  was  10:1. 

lesion  (between  the  lesion  and  the  camera),  and  the  detector 
blur,  both  act  to  reduce  the  lesion  contrast  in  the  image. 

To  illustrate  the  effects  of  attenuation  and  detector  blur,  an 
image  of  a  phantom  containing  2  simulated  lesions  is  shown 
in  Fig.  1.  The  image  shows  2  lesions  of  equal  size  (inner 
diameter  -  1  cm)  and  equal  radioactivity  concentration  (3 
pCi/cc)  but  at  different  depths  in  the  phantom,  and  therefore 
different  distances  from  the  detector.  The  background 
radioactivity  concentration  was  0.3  pCi/cc.  Lesion  1  is  at  a 
distance  of  7.7  cm  from  the  detector  whereas  lesion  2  is  at  a 
distance  of  5  cm.  Even  though  the  lesions  have  equal  size  and 
radioactivity  concentration,  their  image  contrast  is  vastly 
different.  The  image  contrast  for  lesion  1  is  0,26  whereas  that 
for  lesion  2  is  0.42.  This  illustrates  the  fact  that  in  order  to 
quantify  the  radiotracer  concentration,  it  is  necessary  to  make 
corrections  to  the  lesion-to-background  contrast  observed  in 
the  image. 

We  are  developing  a  dual  modality  breast  scanner  with  a 
digital  x-ray  detector  and  a  gamma  ray  detector  mounted  on  a 
single  upright  mammography  unit.  In  the  first  generation  of 
the  scanner,  the  surface  of  the  x-ray  detector  acts  as  the  breast 
support  [4].  Mild  compression  is  applied  and  an  x-ray  and  a 
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gamma  image  are  obtained.  Merging  the  two  images  together 
provides  structural  and  functional  information  in  a  single 
image.  The  second  generation  dual  modality  breast  scanner  is 
designed  to  permit  multiple  views  of  the  breast  with  fixed 
compression.  The  breast  is  held  under  mild  compression  on  a 
breast  support  structure  that  is  independent  of  the  detectors 
and  the  gantry  arm  [5].  The  x-ray  source  and  the  x-ray  and 
gamma  detectors  are  mounted  on  a  motor-driven  gantry  arm. 
The  motorized  gantry  with  the  breast  support  structure  and  a 
small  field  of  view  x-ray  detector  mounted  on  it  is  shown  in 
Fig.  2.  The  mechanical  aspects  of  this  system  have  been 
described  earlier  [5].  Multiple  view  imaging  makes  it  possible 
to  obtain  three-dimensional  lesion  localization  with  this 


Fig.  2.  Photograph  of  the  dual  modality  system  equipped  with  breast 
support  structure  independent  of  the  x-ray  detector.  A  small  field  of  view 
x-ray  detector  is  shown  mounted  on  the  gantry.  However,  the  gamma 
camera  is  not  included  in  the  picture. 

system. 


II.  Stereotactic  Localization 

Lesion  localization  via  stereotactic  imaging  is  a  method 
commonly  used  in  breast  needle  biopsies  [6-8].  The  location 
of  the  lesion  in  the  two  stereo  images  is  used  to  calculate  the 
three-dimensional  location  of  the  lesion  center,  Xq,  yo,  Zo-  The 
expressions  for  Xq,  y^  z©  are  given  in  (l)-(3),  shown  at  the 
bottom  of  the  page. 

X  (fl  -  c)sin  0{(ziX2  +  22-^1 ) **•  (^  ~ g)(sin -z^)-  cos +  X;))) 

(ZjXj -22X,)  +  (fl-c)(sin0(z, +Zj)  +  cos^(Xj -X2))-2(a-c)^sin^cos^ 

(1) 

2  -  ^0^1  ~ (g ~ c)x^ cos^ -(a- c)z^  sin ^  +  (g - c)xi  cos^  (2) 
X, -(a-c)sin^ 

Ly^(x,-(a-c)sm9)  (3) 

"  Lj,,cos^  +  Xj, -(fl-c)sin^ 

The  derivation  of  these  equations  and  definition  of  the  co¬ 
ordinate  system  and  variables  used  are  included  in  the 
appendix. 


III.  Theory 


A.  Background  Concentration 

In  the  gamma  emission  image  of  the  breast  obtained  in  a 
dual  modality  scan,  the  average  number  of  primary  counts 
(i.e.  not  including  scatter)  per  pixel  (Nb)  in  the  background 
region  of  the  image  (i.e.  a  region  of  the  breast  adjacent  to  the 
lesion)  can  be  calculated  by  integrating  the  activity  along  the 
thickness  of  the  breast. 

D 

N^  =  stp^  dz  (4) 

0 

where  s  is  the  detection  efficiency,  t  is  the  acquisition  time,  p 
is  the  pixel  size,  and  Cb(z)  is  the  activity  per  unit  volume  in 
the  background  portion  of  the  breast.  D  is  the  thickness  of  the 
breast  above  the  pixel  and  |it  is  the  linear  attenuation 
coefficient  of  the  breast  tissue  for  140  keV  gamma  rays.  Since 
mild  breast  compression  is  used,  D  is  simply  the  compressed 
breast  thickness.  In  the  case  where  the  background  activity 
concentration  is  uniform  along  the  line  of  integration,  this 
equation  can  easily  be  rearranged  to  calculate  the  background 
concentration. 


(5) 


The  average  number  of  counts  obtained  in  a  background 
pixel  of  the  image  is  the  sum  of  the  scatter  and  the  primary 
counts,  Nb.  These  scatter  counts  must  be  subtracted  from  the 
sum  to  obtain  the  primary  background  counts.  The  scatter-to- 
primary  ratio  (SPR)  is  dependent  on  the  location  within  the 
breast,  increasing  with  decreasing  distance  from  the  chest 
wall  [9].  In  practice,  the  position-dependent  SPR  is  obtained 
by  acquiring  images  in  list  mode.  This  permits  retrospective 
computation  of  position-dependent  pulse  height  spectra  from 
which  the  SPR  can  be  obtained. 


B.  Lesion  Concentration 

There  are  two  primary  factors  that  make  the  contrast  of  the 
lesion  in  the  image  a  poor  measure  of  the  relative  radiotracer 
concentrations  in  the  lesion  and  the  background  breast  tissue. 
These  are  i)  the  attenuation  of  gamma  rays  originating  in  the 
lesion,  caused  by  the  tissue  overlying  the  lesion  and  ii)  the 
blur  due  to  the  detector  point  spread  function  (PSF).  In 
addition,  for  both  lesion  pixels  and  background  pixels,  scatter 
adds  additional  counts  thus  reducing  the  apparent 
concentration  ratio  in  the  image. 

The  lesion  concentration  is  given  by: 

C  ^  total 
VitsA 

where  Nl, total  is  the  total  number  of  detected  counts  from  the 
lesion,  Vl  is  the  lesion  volume,  t  is  the  acquisition  time,  e  is 
the  gamma  camera  efficiency  and  A  is  the  attenuation  of  the 
breast  tissue  overlying  the  lesion.  The  amount  of  attenuation 
is  given  by: 


(7) 
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where  b  is  the  depth  of  the  lesion  in  the  breast. 

The  blurring  due  to  the  detector  PSF  does  not  change  the 
total  number  of  counts  detected  from  the  lesion  but  only 
redistributes  these  counts  over  a  larger  area  thereby  reducing 
the  amplitude  of  the  lesion  peak  in  the  image.  The  projection 
of  the  spherical  lesion  on  a  perfect  detector  (that  is  a  detector 
whose  PSF  is  a  delta  function)  can  be  expressed  as 


=  0/orr  -ol  >1  (8) 

R 

where  R  is  the  lesion  radius,  Xo  and  yo  are  the  x  and  y  co¬ 
ordinates  of  the  lesion  center  and  the  peak  amplitude  NL,max  is 
set  by  the  condition: 


lj^L(x,y)dxdy  =  Ni^^^ 


To  calculate  the  reduction  in  the  amplitude  at  the  lesion 
center  in  the  image  due  to  the  PSF  of  the  detector,  the  two- 
dimensional  projection  of  the  spherical  lesion  on  the  detector 
given  by  (8)  is  convoluted  with  the  two-dimensional  Gaussian 
PSF  of  the  detector,  The  width  of  the  Gaussian  in  both 
dimensions  is  a  function  of  the  distance  of  the  source  from  the 
detector  surface.  The  detector  PSF,  normalized  for  peak 
amplitude  of  unity,  can  be  written  as 

g(x,y)=exp  (10) 

(x,  y)  =  (x,  y)  ®  g(x,  y)  (l  1) 

Convoluting  the  above  two  functions,  we  get  the  blurred 
spatial  distribution  of  the  detected  lesion  counts  in  the  image, 
NL,biur(x,y).  Clearly,  NL,biur(x,y)  is  spread  over  a  larger  area 
and  has  a  reduced  height  relative  to  the  original  lesion 
projection,  NL(x,y).  The  ratio  of  the  peak  amplitude  of 
NL.biur(x,y)  to  that  of  NjL(x,y)  is  defined  as  the  contrast 
reduction  coefficient  (CRC)  [3]. 

=  (12) 

^iL,max 

The  detector  blur  correction  is  applied  by  dividing  the 
amplitude  of  the  lesion  profile  in  the  image,  Ntbiurmax  by  the 
CRC. 

XT  ^Lyblur,m2iX 

(13) 

We  thus  obtain  the  amplitude  of  the  lesion  profile  in  the 
absence  of  detector  blur,  N^.max  and  can  calculate  the  total 
number  of  detected  lesion  counts,  total  from  (8)  and  (9)  and 
finally  the  lesion  concentration,  Cl  from  (6). 


IV.  Phantom  Measurements 

A.  Lesion  Localization  with  x-ray 

A  finable  phantom  of  rectangular  cross  section  was 
designed  to  simulate  a  compressed  breast.  A  tube  was  located 
at  the  phantom  center  to  which  spherical  lesions  were 
attached  using  low  attenuation  rods  of  known  length,  r, 
priented  at  a  known  angle,  cp  with  respect  to  the  phantom 
i!s\ji[face.  A  support  was  designed  to  allow  the  phantom  to  be 
mounted  with  the  tube  on  the  rotation  axis  of  the  gantry.  A 
bearing  allowed  the  gantry  arm  to  rotate  while  the  phantom 
was  held  stationary.  This  setup  enabled  us  to  place  the  lesions 
at  accurately  known  x,  y  and  z  co-ordinates.  Lesions  of 
interior  diameter  0.8  and  1.0  cm  from  Data  Spectrum  Corp., 
Chapel  Hill,  NC  were  used  for  this  study.  The  wall  thickness 
of  the  lesions  was  only  ~  1  mm  and  hence  the  attenuation 
through  the  walls  could  be  ignored  for  the  concentration 
calculation.  With  the  phantom  positioned  independently  of 
the  x-ray  detector,  images  were  acquired  at  several  different 
angles  from  -90  to  +90  degrees.  The  phantom  was  filled  with 
water  and  the  lesions  contained  a  small  amount  of 
gadolinium-DTPA  for  enhanced  x-ray  contrast. 

B.  Concentration  Ratio  Measurement 

For  simplicity,  acquisition  of  the  gamma  images  of  the 
phantom  was  performed  with  the  gamma  camera  positioned 
stationary  on  a  bench  top.  A  schematic  of  the  experimental 
setup  is  shown  in  Fig.  3.  The  lesions  and  background  were 
filled  with  a  solution  of  technetium-99m  with  a  lesion-to- 
background  concentration  ratio  of  10:1  to  simulate  hot  lesions 
in  a  warm  background.  The  absolute  concentration  in  the 
background  was  0.48  pCi/cc  and  that  in  the  lesion  was  4.9 
pCi/cc.  In  order  to  test  our  ability  to  measure  concentration 
ratios  over  a  range  of  lesion  depths,  the  central  tube  of  the 
phantom  was  attached  to  a  rotation  stage  from  Velmex,  Inc., 


e  ^"7"^ 


Fig.  3.  Experimental  setup  used  for  the  phantom  measurements.  A 
rectangular  box  phantom  was  used.  The  central  tube  with  lesions  attached 
by  low  attenuation  rods  was  rotated  independent  of  the  outer  box.  Based 
on  the  known  values  of  r  and  (p,  the  three-dimensional  location  of  each 
lesion  for  each  0  value  could  be  calculated. 


^  Paper  #M6-237 


Bloomfield,  NY  (website:  http://www.velmex.com),  and 
rotated  independently  of  the  fixed  phantom  housing  over  a 
range  of  360  degrees.  At  each  angle,  0,  the  co-ordinates  Xo, 
Yo,  Zo  of  the  lesion  centers  were  known  based  on  the  known 
rod  lengths  and  angular  spacing  between  lesions.  Images  were 
acquired  every  15  degrees  over  a  range  of  360  degrees. 
Rotating  the  tube  with  respect  to  the  phantom  housing  and 
detector  surface  places  the  lesions  at  varying  depth  and 
lesion-to-collimator  distance,  thereby  permitting  the 
correction  model  to  be  tested  over  a  range  of  lesion  depths 
with  a  single  setup. 

V.  Detector  Characterization 

The  gamma  camera  used  for  the  studies  described  here 
consists  of  a  4  X  4  array  of  R7600-00-C8  position  sensitive 
photomultiplier  tubes  (PSPMTs)  from  Hamamatsu, 
Bridgewater,  New  Jersey.  The  PSPMTs  are  optically  coupled 
via  acrylic  light  guides  to  a  30  x  30  array  of  Nal(Tl)  crystals 
from  Saint-Gobain,  Newbury,  Ohio.  Each  crystal  is  3  mm  x  3 
mm  X  6  mm  (thick),  and  the  crystal  center-to-center  spacing  is 
3.3  mm.  The  active  field  of  view  of  the  gamma  camera  is  ~  10 
cm  X  10  cm  and  the  'energy  resolution  is  15%  at  140  keV.  A 
lead  collimator  with  1.58  mm  hexagonal  apertures  and  height 
of  21  mm  was  used.  The  septal  thickness  was  0.267  mm.  An 
energy  window  of  -7%/+22%  centered  around  the  140  keV 
photopeak  was  used.  The  camera  sensitivity  with  this  energy 
window  is  220  cpm/pCi  (efficiency  of  1.0  x  lO"*).  Further 
details  on  the  detector  design  and  performance  are  published 
elsewhere  [10], 

The  full  width  at  half  maximum  (FWHM)  of  the  two- 
dimensional  gamma  camera  PSF  was  obtained  by  measuring 
the  image  of  a  *’Co  spot  marker  from  Isotope  Products, 
Valencia,  CA  at  a  variety  of  source-to-collimator  distances. 
For  each  source-to-collimator  distance,  the  spot  marker, 
which  has  an  active  diameter  of  ~1.5  mm,  was  imaged  at 
several  locations  on  the  detector  surface.  These  images  were 
shifted  and  summed  to  obtain  a  source  image  profile  averaged 
over  the  camera  surface.  The  profile  of  the  point  source, 


which  was  assumed  to  be  a  uniform  disc  of  diameter  1.5  mm, 
was  deconvoluted  from  the  averaged  profile.  The  result  was 
then  fitted  with  a  two-dimensional  Gaussian  function.  The 
FWHM  of  the  fitted  function  in  the  x-  and  y-dimensions  are 
plotted  as  a  function  of  the  source  to  detector  (source  to 
crystal)  distance  in  Fig.  4. 

VI.  Data  Analysis 

A.  Lesion  Localization 

An  analysis  program  written  in  IDL  (Interactive  Data 
Language  from  Research  Systems,  Inc.,  Boulder,  CO, 
http://www.rsinc.com)  was  used  to  locate  the  lesion  center  in 
the  stereo  images.  Using  (1),  (2)  and  (3)  the  co-ordinates  of 
the  lesion  center  were  calculated.  The  z  value  of  the  lesion 
was  then  used  to  calculate  the  lesion  and  background 
concentration  from  the  image  counts. 

B,  Concentration  Ratio  Measurement 

1)  Background  Concentration:  To  get  the  number  of 
background  counts,  a  profile  through  the  gamma  image  is 
taken  which  does  not  contain  any  lesions.  For  each  lesion, 
profiles  are  extracted  from  rows  near  that  lesion.  From  these 
rows,  the  average  number  of  background  counts  is  obtained 
using  pixels  taken  from  columns  containing  the  lesion. 

‘ihe  scatter  contribution  to  this  background  pixel  count  is 
subtracted  to  obtain  the  primary  counts  in  the  background. 
For  the  phantom  used  here,  the  scatter  fraction  is 
approximately  constant  and  uniform  throughout  the  region  in 
which  the  lesions  are  located.  The  technique  for  obtaining  the 
scatter  fraction  for  a  box  phantom  such  as  that  used  in  these 
experiments  has  been  described  earlier  [5].  The  scatter 
fraction  for  the  energy  window  -7%/+22%  was  measured  to 
be  0.2.  The  background  concentration  is  then  determined 
from  the  primary  background  counts  as  described  in  (5). 

2)  Lesion  Concentration:  To  calculate  the  lesion 
concentration,  a  region  of  interest  centered  on  the  lesion  is 
first  extracted.  The  number  of  counts  in  a  pixel  imaging  the 
lesion  is  the  sum  of  lesion  counts  and  counts  from  the 
background  volume  above  and  below  the  lesion.  The  average 
number  of  background  counts  must  be  subtracted  from  the 
total  number  of  lesion  pixel  counts  to  obtain  the  lesion-only 
counts.  In  order  to  obtain  a  good  estimate  of  the  amplitude  of 
the  lesion  profile  in  the  presence  of  statistical  fluctuations,  a 
two-dimensional  surface  is  first  fitted  to  the  background 
subtracted  lesion  profile.  This  lesion  profile  amplitude  must 
be  divided  by  the  CRC  to  compensate  for  the  detector  blur. 

To  obtain  the  CRC  the  following  procedure  is  used.  The 
known  lesion-to-collimator  distance  is  used  to  determine  the 
X-  and  y-widths  of  the  Gaussian  detector  PSF  (Fig.  4).  The 
PSF  is  convoluted  with  the  theoretical  lesion  projection  as  in 
(11)  and  the  theoretical  CRC  is  obtained.  The  known  lesion 
radius  was  used  for  the  value  of  R  in  (8).  The  peak  height 
obtained  after  dividing  the  fitted  lesion  amplitude  by  the  CRC 
is  then  attributed  to  the  maximum  of  the  unblurred  lesion 
projection  function,  NL(x,y),  and  the  area  under  it  is 
calculated  to  obtain  the  total  counts  coming  from  the  lesion. 
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These  counts  are  then  divided  by  the  attenuation  of  the 
overlying  tissue.  For  the  case  of  the  phantom  study  reported 
here,  the  overlying  tissue  is  water  and  the  phantom  wall.  The 
attenuation  of  the  gamma  rays  from  the  lesion  is  given  by: 

A^^ll  (14) 

where  n«,  is  the  linear  attenuation  coefficient  of  water  for  140 
keV  gamma  rays  and  A^an  is  the  attenuation  of  the  acrylic 
wall,  on  the  camera  side  of  the  phantom. 

VII.  Results 

A.  Lesion  Localization 

Table  1  shows  the  results  obtained  for  the  stereotactic 
localization  using  the  x-ray  imaging  component  of  the  dual 
modality  system.  These  were  obtained  with  stereo  views  at 
angles  +15  degrees  and  -15  degrees.  Columns  2  and  3  give 
the  known  location  of  the  lesion  center  in  the  x  and  z 
dimensions  and  columns  4  and  5  give  their  location  as 
calculated  from  the  stereo  views  and  equations.  The  last 
column  gives  the  error,  obtained  by  subtracting  the  calculated 
values  from  the  known  values  in  millimeters.  Similar  results 
were  obtained  using  other  stereo  angles  up  to  ±60®. 


Table  I 

Measured  and  Calculated  Lesion  Locations 


Lesion  # 

Known 
xo  (cm) 

20  (cm) 

Calculated 

xo(cm) 

20(cm) 

Difference 
X  (mm) 

2  (mm) 

1 

nofn 

■ilSM 

|Q| 

Bhi!I 

0.05 

1^9 

B.  Concentration  Ratio  Measurement 

Fig.  5  shows  the  ratios  between  the  concentrations  in  two 
of  the  simulated  lesions  and  the  background,  plotted  as  a 
function  of  the  rotation  angle  of  the  central  tube.  Both  lesions 
had  an  inner  diameter  of  1.0  cm  and  a  lesion  to  background 
ratio  of  10:1.  Plots  are  shown  of  concentration  ratios 
calculated  both  with  and  without  corrections  for  detector  blur, 
attenuation,  and  scatter.  Corrected  ratios  were  calculated 
using  the  method  described  above.  Uncorrected  concentration 
ratios  were  obtained  for  each  view  (all  0  values)  as  follows. 
ROIs  were  drawn  containing  the  lesion  image,  and  in  a  region 
nearby  containing  only  background.  The  average  background 
count  per  pixel  was  subtracted  from  each  pixel  in  the  lesion 
ROI.  A  smooth  two-dimensional  function  was  fitted  to  the 
resulting  lesion  profile,  and  its  area  was  calculated.  This  area, 
divided  by  the  known  lesion  volume,  was  taken  as  the  lesion 
concentration.  The  background  concentration  was  taken  to  be 
the  total  number  of  counts  in  the  background  ROI  divided  by 
the  product  of  the  ROI  area  and  the  phantom  thickness.  Fig.  5 
shows  that  with  corrections  made,  the  concentration  ratios  are 
closer  to  the  true  concentration  ratio. 

Plot  6  shows,  for  every  0,  the  image  contrast  obtained 
using  the  measured  lesion  and  background  pixel  counts. 


without  any  corrections  made  for  attenuation  or  resolution 
degradation  (squares).  This  illustrates  that  estimation  of 
concentration  ratios  based  simply  on  image  contrast  can  be 
subject  to  large  errors. 


plotted  as  a  function  of  the  rotation  angle  for  two  lesions  with  different  r 
(r  =  1.4  cm  in  the  upper  graph  and  r  =  2.7  cm  in  the  lower  graph).  The 
difference  in  the  depth  variation  (2.8  cm  versus  5.4  cm)  is  reflected  in 
the  difference  in  die  shapes  of  the  uncorrected  curves. 


concentration  ratio  is  shown  in  the  upper  graph  of  Fig.  5. 
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VIIL  Discussion  and  Conclusions 

The  results  obtained  for  the  stereotactic  lesion  localization 
demonstrate  accuracy  to  within  1-2  mm.  This  is  accurate 
enough  to  be  able  to  use  the  information  to  make  corrections 
for  the  attenuation  and  detector  blur  and  to  significantly 
improve  our  estimate  of  the  true  concentration  ratio  in  the 
tumor  and  the  background,  as  is  demonstrated  in  Fig.  5. 
However,  for  this  phantom  study,  our  calculations  used  the 
known  lesion  radius  and  volume.  In  the  clinical  situation, 
these  parameters  are  unknown,  and  must  be  estimated  from 
the  images.  Hence,  methods  need  to  be  developed  to  extract 
this  information  from  the  high-resolution  x-ray  image. 
Similarly,  the  scatter  fraction  in  the  breast  is  a  parameter  that 
must  be  measured.  The  scatter  fraction  is  a  function  of  the 
location  in  the  breast,  and  increases  with  decreasing  distance 
to  the  chest  wall  [9].  Position-dependent  energy  spectra  are 
routinely  obtained  during  a  dual  modality  scan,  and  can  be 
used  to  obtain  the  position  dependent  scatter  fraction. 

In  this  paper  we  have  described  lesion  localization  using  x- 
ray  imaging.  In  principle,  the  three-dimensional  lesion 
location  could  also  be  obtained  using  gamma  ray  stereo 
imaging  [11].  However,  gamma  image  acquisition  times  are 
significantly  longer  than  those  of  x-ray  images,  and  the 
accuracy  of  localization  is  lower  due  to  the  lower  resolution 
of  the  gamma  detectors.  Also,  it  is  desirable  to  have  two  x-ray 
views  in  order  to  better  estimate  the  lesion  shape  and  volume. 
We  have  therefore  elected  to  use  stereotactic  x-ray  imaging 
for  lesion  localization  rather  than  gamma  imaging.  However, 
in  situations  where  the  lesion  is  poorly  visualized 
radiographically  (e.g.  radiodense  breasts),  but  is  well 
visualized  in  the  gamma  images,  stereotactic  gamma 
localization  would  be  the  approach  of  choice. 

Our  phantom  studies  suggest  that  with  appropriate 
corrections  for  attenuation,  detector  blur  and  scatter,  based  on 
measured  lesion  location,  we  can  more  accurately  predict 
radioactivity  concentration  ratios.  Stereotactic  x-ray  imaging 
combined  with  gamma  ray  imaging  is  a  potentially  powerful 
method  to  quantify  radiotracer  uptake  in  breast  tumors. 

IX.  Appendix 

Following  is  the  derivation  of  equations  (1)  -  (3),  which 
calculate  the  co-ordinates  of  the  lesion  center  based  on  two 
stereo  views.  The  co-ordinate  system  is  shown  in  Fig.  7.  We 
define  the  origin  of  the  co-ordinate  system  to  be  the  point  on 
the  rotation  axis  that  is  intersected  by  the  line  between  the  x- 
ray  focal  spot  and  the  center  point  of  the  chest  wall  edge  of 
the  x-ray  detector.  The  x-axis  is  defined  to  lie  along  the 
medial-to-lateral  direction  of  the  patient,  with  positive  x 
pointing  towards  the  right  of  the  patient.  The  y-axis  is  along 
the  rotation  axis  of  the  gantry  with  the  positive  direction 
going  towards  the  gantry  and  away  from  the  patient.  The  z- 
axis  is  defined  to  be  along  the  detector  to  source  direction 
when  the  gantry  arm  is  positioned  vertically,  with  positive  z 
going  towards  the  source.  Let  a  be  the  source  to  detector 
distance  and  c  be  the  distance  between  the  origin  and  the 
center  point  of  the  chest  wall  edge  of  the  detector  surface. 


source,  S,  when  the  gantry  arm  is  vertical  (0  =  0)  and  (xj,  yd, 
Zd)  are  the  co-ordinates  of  the  detector  chest  wall  edge  center, 
D,  when  0  =  0.  O  (x,  y,  z)  is  the  projection  of  the  center  of  the 
lesion  with  spatial  co-ordinates  (Xo,  yo,  Zo)  on  the  detector 
plane. 

If  the  source  and  detector  are  rotated  through  an  angle  0i , 
the  new  source  and  detector  center  co-ordinates  are  S'(Xsi,  ysi, 
Zsi)  and  D^Xdi,  ydi,  Zdi)  respectively,  and  the  image  of  the 
lesion  center  on  the  detector  is  located  at  O^Xi,  yi,  Zi).  We 
denote  the  perpendicular  distance  in  the  images  between  the 
lesion  center  and  the  projection  of  the  rotation  axis  onto  the 
detector  by  Lxi,  and  denote  the  perpendicular  distance 
between  the  lesion  center  and  the  chest  wall  edge  of  the 
detector  by  Lyi.  From  geometry  we  get  the  following 
relationships, 

^.1  =(«-c)sin/9,  ;  =0  ;  =(a-c)cos(9, 

XJ^=-csme^  ;  7^1=0  ;  =-ccos^, 

x^  =  4,  cos  (9,  +  ;  7,  =  4,  ;  -  4,  sin  9^ 


The  equation  of  line  S'O'  is  given  by, 

~  _  -To  ~  _  ^0  ~ 

~  yi  ~  j'ji  ^1  “  ■2,1 
x^-(a-c)sm9^  z„ -(a-c)cos^i 
X, -(a-c)sm0i  z, -(a-c)cos^, 
Let  R  =  (a  -  c)sm6^ 
andT  =  {a-c)cos6^ 

then; 

-R  _  Zq -t 

x^-R  z^-T 

.^AXo-R){z,-T) 

JCj  -  7? 


(15) 

(16) 


(17) 

(18) 


P^er  #M6~237 


For  a  negative  rotation  through  an  angle  02,  the  new  source 
and  detector  center  co-ordinates  are  S"(Xs2,  ys2,  Zs2)  and 
D"(Xd2,  Ydi,  Zd2)  respectively.  O"  (X2,  yi,  Z2)  is  the  projection 
of  the  lesion  center  on  the  detector  plane  in  this  configuration, 
We  now  have  the  following  equations, 

x,2=(a-c)sin02  ;  7.2=0  ; 

xj2=-csm02  ;  yd2=0  ; 

x2=L^2<^o^^2+Xd2  ;  y2=h2  J 

I'he  equation  of  line  S”0”  is  given  by, 

-^.2  _  7o  -7.2 

X. -X 


z,2  =(<2-c)C0S^2 

Zd2  =-C  008^2 


z.  -z 
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"^2-^.2  72-7.2 


^2-2.2 


(19) 


-ja-c) sin ^2  z^-{a-c)cos6^ 


^2 -(a- c)  sin  6*2  Z2 -(a-c)cos6»2 
Lets  =  (a-c)sin02 
andU  =  (a  -  c)  cos  ^2 
z-U 


(20) 


;c2-5 


z^-U 


.  _fe-^)(Z2-^) 
^  x^-S 


(21) 


(22) 


z„  =■ 


■  +  T 


(24) 


We  can  solve  equations  (18)  and  (22)  simultaneously  for  Xo 
and  Z(,.  We  get; 

^  _  R(.2,  -  T){x,  -S)-  S(z,  -  U)(x,  -R)  +  (T-  UXx,  -  R)(x,  -  5') 

°  -  (zj  -  U)(x,  -R)  +  (z,  -  T)ix,  -  S) 

(23) 

_(x„-R)(z,-T) 

(x^-R) 

If  the  two  stereo  views  are  chosen  to  be  symmetric  about  the  z 
axis  then  0|  =  -02  =  0  and 

COS  =  COS  6*2  =  cos  0  and  sin(9,  =-sin^2  =sin^ 

We  now  substitute  for  R,  S,  T,  U  to  get 

^  (a-c)5ing((z,j:,  +z,jC|)  +  (a-eXsing(2|  -z,)-cosg(x, 

(z, Jtj  -  ZjX, )  +  (o  -  c)(sin  0{z,  +  z, ) + cos  0(x,  -xj)-  2(a  -  c)‘  sin  0  COS  ^ 

(25) 

^  cos^-(fl-c)Z|  s^n^  +  (fl-c)x^  cos^ 

X, -(a-c)sin^ 

.  t 

To  solve  for  yoi 

yi-ysi 

x,-(a-c)sm0  K 

4,cos(9+x^,-(a-c)sin(9 
^y,(x„-(a-c)sin&) 


>'»  = 


L^j  cos0  +  Xj,  -(a-c)sm0 


(26) 

(27) 

(28) 

(29) 


Equations  (25),  (26)  and  (29)  are  equations  (l)-(3)  in  Section 

11. 
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